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ABSTRACT 

Starting from the SFR per unit physical area, determined for DLAs using the C II* method, we obtain 
the SFR per unit comoving volume at z « 3. Pure warm neutral medium (WNM) solutions are ruled 
out since they generate more bolometric background radiation than observed, but the CNM-dominated 
two-phase solutions are consistent with the backgrounds. We find the SFR per unit comoving volume 
for DLAs agrees with that for the Lyman Break Galaxies (LBGs). Though the mass of produced stars 
indicated by the SFRs is consistent with the current densities of known stellar populations, the mass of 
metals produced by z=2.5 is 30 times larger than detected in absorption in DLAs. Of the three possible 
solutions to this "missing metals" problem, the most likely appears to be that star formation occurs in 
compact bulge regions. We search for evidence of feedback and find no correlations between the SFR per 
unit area and N(H I), but possible correlations between SFR per unit area and low-ion velocity width 
and SFR per unit area and metal abundance. We show that (a) the correlation between cooling rate and 
dust-to-gas ratio is positive evidence for grain photoelectric heating, (b) the CMB does not significantly 
populate the C II excited fine-structure states, and (c) the ratio of CII* to resonance-line optical depths 
is a sensitive probe of the multi-phase structure of the DLA gas. We address recent arguments that 
DLAs are comprised only of WNM gas, and show them to be inconclusive. Despite the rough agreement 
between SFR per unit comoving volume for DLAs and LBGs, current evidence indicates these are distinct 
populations. 

Subject headings: cosmology — galaxies: evolution — galaxies: quasars — absorption lines 



1. INTRODUCTION 

This is the second of two papers describing a new method 
for obtaining star formation rates (SFRs) in damped Lya 
systems (DLAs) . In Paper I (Wolfe, Prochaska, & Gawiser 
2003) we showed how measurements of C II* 1335.7 ab- 
sorption lines in DLAs allow one to infer the cooling rate 
per H atom of the neutral gas. Since we assume steady- 
state conditions, this equals the heating rate per H atom, 
which we use to infer the SFR per unit area, "0* • We do this 
by assuming gas in DLAs to be heated by the same mech- 
anism responsible for heating the ISM in the Milky Way, 
the grain photoelectric effect (Bakes & Tielens 1994; Wein- 
gartner & Draine 2001a). In that case the heating rate is 
proportional to the product of the dust-to-gas ratio, k, the 
photoelectric heating efficiency, e, and the mean intensity 
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vatory is a joint facility of the University of California and the Cali- 
fornia Institute of Technology. 

^Current address: Yale Astronomy Department, P. O. Box 
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of FUV radiation, Jj^; the latter is proportional to for 
sources in a plane parallel layer. Specifically, in Paper I we 
modeled DLAs as uniform gaseous disks with radius, R, 
and scale-height h, in which the sources of FUV radiation 
were uniformly distributed. We also showed how k can be 
deduced from the [Fe/Si] and [Si/H] abundance ratios (re- 
cah [X/Y]EElogio(X/Y) -logio(X/Y)Q) for the following 
assumptions about grain composition: grains were either 
carbonaceous as in the Galaxy (the "Gal" model) or Sil- 
icates as in the SMC (the "SMC" model). We inferred k 
by assuming the number of depleted C or Si atoms per 
depleted Fe atom to be the same in DLAs as in the ISM. 
Furthermore, we considered depletion ratios ranging from 
a minimal "nucleosynthetic ceiling" in which the intrin- 
sic ratio, [Fe/Si]i„t=— 0.2, to a maximal depletion ratio, 
[Fe/Si]„,t=0.0 (see Prochaska & Wolfe 2002; hereafter re- 
ferred to as PW02). 

In Paper I we solved the transfer equation for and 
then calculated the thermal equilibrium of gas subjected to 
cosmic-ray and X-ray heating in addition to grain photo- 
electric heating. The gas was assumed to cool in the usual 
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way; i.e., by emission of fine-structure, metastable, and 
Lya lines as well as grain recombination radiation. We 
found that gas can reside in two thermally stable states; a 
cold neutral medium (CNM) and a warm neutral medium 
(WNM) (sec Wolfirc ct al. 1995, hereafter W95). Typ- 
ically, the densities and temperatures of the CNM and 
WNM are 10 cm'^and 150 K, and 0.2 cm'^ and 8000 K 
respectively. We further assumed the CNM and WNM to 
be in pressure equilibrium at pressure Peq — {PminPmax)^^'^ 
where Pmin and Pmax arc the minimum and maximum 
pressures of the pressure versus density curve. We con- 
sidered a CNM model in which the line-of-sight to the 
background QSO encounters comparable column densi- 
ties of gas in the CNM and the WNM. We also consid- 
ered a WNM model in which the line-of-sight encoun- 
tered only WNM gas at pressure equal to Peg- Combining 
the measured heating rates with those predicted at the 
thermally stable gas densities, ucnm and hwnm, wc ob- 
tained unique values for for each DLA; one value for the 
CNM solution and the other for the WNM solution. We 
then averaged ?/'*, for two redshift bins centered at z=2.15 
and z=3.70, to derive the average SFR per unit physical 
area, < ^p*{z) >. The WNM models result in significantly 
higher SFRs than the CNM models since the measured 
[C II] 158 /im cooling rate per H atom, Ic, is a small fraction 
of the total cooling rate in the WNM, whereas Ic equals 
the total cooling rate in the CNM (see Paper I). 

This paper starts by considering quantities with cos- 
mological significance. Specifically, in § 2 we combine 
< 4^t:{z) > with the incidence of DLAs per unit absorp- 
tion distance interval, dN/dX (Bahcall & Peebles 1969), 
to derive the SFR per unit comoving volume, p,:{z). We 
then derive the bolometric background intensity, Iebl, for 
all model combinations. We show that the WNM models 
produce more background radiation than observed in every 
case, and as a result are ruled out. By contrast the CNM 
models are consistent with the observed values of Iebl- 
We compute a consensus model, which is an average over 
all the CNM models. We show that the resulting p'*(z) arc 
comparable to inferred for the Lyman Break Galax- 

ies (Steidel et al. 1999; hereafter referred to as LBGs). 
In § 3 we consider implications of these results. We com- 
pute the mass of stars and the mass of metals produced by 
the star formation history, of the consensus model. 

While the mass of stars is consistent with masses of current 
stellar populations, the mass of metals produced by z=2.5, 
is more than 30 times the mass of metals inferred for DLAs 
at the same redshift. We discuss possible solutions to this 
dilemma including a "bulge" model in which star forma- 
tion is confined to a compact region located at the cen- 
ter of the extensive region creating C II* absorption. We 
consider independent evidence for (a) star formation and 
(b) the deposition of stellar energy into the absorbing gas, 
i.e., feedback. At this point the reader not interested in 
the physics of interstellar gas can turn to the final section, 
§ 6. Having discussed various implications of our mod- 
els wc proceed to test their self-consistency in § 4 where 
three tests are carried out. First, we find a statistically 
significant correlation between the [C II] 158 /xm cooling 
rate per atom, Ic, and k, which is strong evidence in favor 
of grain photoelectric heating. Second, we show that the 
spontaneous energy emission rate, Ic, reflects the cooling 



rate of the gas instead of the excitation level caused by 
CMB radiation. Third, we examine the ratio of C II* to 
resonance-line optical depths to look for evidence of shifts 
in gas phase. In § 5 we discuss arguments made by other 
authors against the presence of CNM gas in DLAs. A care- 
ful reassessment of these arguments shows that they do no 
rule out the presence of CNM gas in DLAs. A summary 
and concluding remarks are given in § 6. 

Unless stated otherwise we adopt an Einstein-deSitter 
cosmology in which f7M=l, nA=0, and h=0.5 to facilitate 
comparison with published results. 

2. COSMOLOGICAL QUANTITIES 

We now turn to quantities with cosmological signifi- 
cance. We compute the SFR per unit comoving volume, 

p*. The redshift dependence of implies a star forma- 
tion history throughout spacetime that gives rise to back- 
ground radiation. We calculate the bolometric intensity of 
this background radiation for the CNM and WNM mod- 
els and compare the results with observations. We then 
construct a consensus CNM model for /9'*(z), which is con- 
sistent with measurements of the background radiation. 

2.1. The SFR per Unit Comoving Volum,e 
The SFR per unit comoving volume for DLAs is given 



by 



p*{z) = < ip*{z) >Anco{z) 



(1) 



where < > is the average SFR per unit physical area 

at redshift z, and A and the average physical cross- 

sectional area and comoving density of the DLAs. While 
neither A nor Uco has been determined from observations, 
their product is known from the incidence of DLAs per 
unit absorption distance interval, dA/" /dX (e.g. Storrie- 
Lombardi & Wolfe 2000). We find 

dM/dX ^ ApUcoiX) , (2) 

where Ap is the average projection of A on the plane of 
the sky, and X{z), the absorption distance (Bennett et al. 
2003), is given by 



dX 
dz 

As a result 



[(1 + z)^{l + Qmz) - z{z + 2)QaV/^ 



p;(z) = < ^,{z) >{A/Ap)dU/dX 



(3) 

(4) 

We computed by assuming the DLAs to be plane- 
parallel layers; i.e., A/Ap=2, and by choosing an Einstein- 
deSitter cosmology (nM=l,f^A=0, /i=0.5). Although this 
model is ruled out by observations (e.g. Bahcall et al. 
1999), it is the model used in most published determina- 
tions of p« (e.g., Steidel et al. 1999; Lanzetta et al. 2002), 
and for this reason we selected it for purposes of compar- 
ison. We chose < V'*('Z) > and its associated errors from 
Table 3 in Paper I and calculated dN jdX at the mean z 
of the redshift bins from the expression 

dAr/dX = (J7o/c)0.055x(l-f z)°-6i , (5) 

found by Storrie-Lombardi & Wolfe (2000) for the Einstein- 
deSitter cosmology. Errors in were computed by prop- 
agating errors in ^* and in dM jdX. 
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Fig. 1. — SFRs per unit comoving volume for "Gal" dust model 
shown as magenta data points for minimal depletion and green data 
points for maximal depletion. Results for CNM model in Figure la 
and for WNM model in Figure lb. Red dotted data points depict 
p* inferred from galaxy luminosities (cf. Steidel et al. 1999; Lilly 
et al. 1996; and Barger et al. 2000). Solid blue curves are fits to 
high-z DLA data that also agree with galaxy data at z <1. Dashed 
blue curves are fits to medians of DLA SFRs and are extrapolated 
to low 2 assuming < tpt,{z) > is constant and combining equations 
(5) and (4). 



The results for the "Gal" dust models are shown in Fig- 
ure 1 as magenta data points for minimal depletion 
and green data points for maximal depletion. As expected, 
the SFRs per unit comoving volume for the WNM mod- 
els (Figure lb) are at least 10 times higher than for the 
corresponding CNM models (Figure la). Furthermore, for 
every model, minimal depletion gives rise to higher than 
the maximal depletion. This is because for a given heating 
rate, is inversely proportional to the dust-to-gas ratio, 
K, and K is lower for minimal depletion than for maximal 
depletion. Figure 1 also reveals no evidence for statisti- 
cally significant redshift evolution of determined by the 
C II* technique. This is in accord with determinations of 

from luminosities measured for flux-limited samples of 
galaxies, shown as red data points. For the galaxy sam- 
ple the SFRs in the two highest redshift bins are based on 
Lyman-break galaxies that are luminous at rest-frame UV 
wavelengths (Steidel et al. 1999), while the four lowest 
redshift points are based on galaxies with lower luminosi- 
ties (Lilly et al. 1996). The bin at 2: = 2 is based on 
far-infrared (FIR) luminous galaxies detected by SCUBA 
(Holland et al. 1999). The redshifts for these objects were 
determined from a still uncertain radio-FIR photometric 
redshift indicator (Barger et al. 2000). Interestingly, the 
magnitude of the comoving SFRs deduced by the C II* and 
galaxy luminosity techniques are not very different in the 
redshift interval where they overlap; i.e., 2:«[2,4.5]. For 
the CNM model the difference is less than a factor of 2, 
while for the WNM model the difference is about a factor 
of 10. We discuss possible implications of this agreement 
in § 2.3 and § 6. 

To test the generality of these conclusions, the calcu- 
lations were repeated for the "SMC" dust models. The 
results, shown in Figure 2, reveal the same patterns as 
found for "Gal" dust. The principal difference is the sys- 
tematically higher values of p* predicted by corresponding 
"SMC" models. This is because the photoelectric heating 
efficiency of silicate grains is lower than that of carbona- 
ceous grains (see Figure 15 in Weingartner & Draine 



Fig. 2. — Same as Figure 1, except "SMC " dust model is assumed. 



2001a), and as a result higher SFRs are required to achieve 
a given heating rate in the case of silicate (i.e., "SMC") 
grains. Comparison between Figures 1 and 2 also re- 
veals smaller differences between yo* derived for minimal 
and maximal depletion in the case of "SMC" dust than 
for "Gal" dust. The phenomenon is present in the CNM 
model but not the WNM model. In the case of maximal 
depletion, [C/H]gos is larger by 0.2 dex than for minimal 
depletion. But an increase in [C/HJ^qs causes a decrease 
in ncNM, which can increase significantly if ncNM in- 
tersects the rising Icrin) curves caused by the transition 
from WNM to CNM temperatures (see Figure 5 in Paper 
I). In the case of "SMC" dust this effect compensates for 
the decrease in -0, caused by the increase in k discussed 
above, because the heating rate Td increases more rapidly 
with density than for "Gal" dust (Weingartner & Draine 
2001a). Consequently, the net difference in -0*, hence p^,, 
is smaller for "SMC" than for "Gal" dust. 

2.2. Bolometric Background Intensity 

While the general trends in the p* versus z plane ap- 
pear to be insensitive to our choice of model assumptions, 
large systematic uncertainties in the comoving SFRs re- 
main. Though uncertainties in dust composition and de- 
pletion level are contributing factors, the largest source of 
error stems from uncertainties in the thermal phase: Does 
C II* absorption arise from gas in the CNM or WNM 
phase? To address this question we compute the extra- 
galactic bolometric background intensity, Iebl- Because 
Iebl is generated by a given star formation history, and 
since the star formation histories indicated by the CNM 
and WNM models are very different, measurements of the 
background intensity may be able to discriminate between 
them. 

The bolometric extragalactic background intensity gen- 
erated by a given star formation history, p*{z), is as fol- 
lows: 

iEs. = ^l^^J^Mt~mt')dt' , (6) 

where L{t') is the bolometric luminosity per unit mass as 
a function of age t' of a stellar population with a specified 
IMF (Madau & Pozzetti 2000), ti? is the formation epoch 
of the stellar population, and is the current age of the 
universe. We tested the models by computing backgrounds 
generated by analytic fits to the C II* comoving SFRs. The 
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fits are shown as smooth curves in Figures 1 and 2, and 
the backgrounds they generate as corresponding curves in 
Figure 3. 

First, we explored the hypothesis that high-z DLAs evolve 
into normal low-z galaxies. Evolution into normal galaxies 
is consistent with the recent identification of DLAs at z < 
1 with the local galaxy population (Zwaan et al. 2002; 
Turnshek et al. 2002; Rosenberg & Schneider 2003). In 
that case the star formation histories are constrained to 
pass through the comoving SFRs measured for galaxies at 
z < 1 and for DLAs at higher redshifts. The results are 
shown as solid blue curves for "Gal" dust in Figure 1 and 
"SMC" dust in Figure 2. For the CNM model (Figures la 
and 2a) the curves resemble the star formation histories 
inferred from galaxy luminosities in the redshift interval z 
= [0,5]. For the WNM model (Figures lb and 2b) the fits 
greatly exceed the comoving SFRs inferred for the LBGs 
at z > 2, though they are in good agreement with the 
galaxy data at 2 < 1. Here the DLAs could represent 
a population of objects undetected in emission at high z 
that evolve into normal galaxies at low redshifts. Second, 
we considered an hypothesis in which the star formation 
histories of DLAs are dictated solely by their redshift evo- 
lution, without regard to the galaxy data. In this case we 
combined the expression for p'*(z) in equation (4) with the 
expression for dJ\f/dX in equation (5). Although we find 
<ip^,{z) > is independent of redshift in the redshift interval 
z = [1.6,4.5], at lower redshifts there are no measurements 
of C II*, and as a result <ijji.{z)> is unknown. Low red- 
shifts are crucial because that is where most of the back- 
ground radiation arises. For simphcity we let <■)/'* (z) > 
equal a constant evaluated by averaging over all the ip^ 
inferred from C II* absorption. The results are shown as 
dashed curves in Figures 1 and 2. 

The resulting backgrounds are shown in Figure 3 where 
we plot the predicted bolometric intensity from DLAs with 
z > Zmin versus Zmin- The significance of this quantity, 
lEBL{z>Zmin), is that it rcvcals the contribution to the 
measured background, Iebl (i-e, Iebl{z>0)), made by 
DLAs in the redshift range for which (z) has been deter- 
mined; i.e., z=[1.6,4.5]. The backgrounds were obtained 
with an Einstein-deSitter cosmology assuming h = 0.5, by 
adopting the same IMF used to relate ■0* to mean intensity 
(see Paper I), using Bruzual and Chariot's (1993) popula- 
tion synthesis libraries (see Madau & Pozzetti 2000), and 
by assuming a formation redshift, zp = 5. For compar- 
ison, the two horizontal straight lines depict upper and 
lower limits on Iebl set by measurements between 0.15 
/im and 1000 /im (Hauser & Dwck 2001). This wavelength 
range is relevant since it brackets the background spectra 
predicted for most models of DLAs (e.g. Pei et al. 1999). 
According to Dwek (2002) the upper limits, which are cru- 
cial here, are conservative and should be regarded as 95 % 
confidence limits. It is important to emphasize that while 
p*{z) is cosmology dependent, the backgrounds computed 
from the C II* technique (i.e. by combining equations (4) 
and (6)) are independent of the adopted cosmology and 
Hubble constant. Therefore, consistency between theory 
and observation is independent of the cosmology one assumes. 

The background intensities predicted for the WNM mod- 
els are too high. Figures 3b and 3c depict backgrounds 
generated by the "Gal" and "SMC" dust models 
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Fig. 3. — Smooth curves depict bolometric background intensity 
due to DLAs with z > z^in- Horizontal lines are empirical 95 % 
confidence upper and lower limits. Figure 3a shows solutions for 
CNM model. Dashed curves labeled "G" and "S" are generated by 
"Gal" and "SMC" star formation histories shown as dashed curves 
in Figures la and 2a. Solid curves so labeled are generated by "Gal" 
and "SMC" star formation histories shown as solid curves in Figure 
la and 2a. Figure 3b shows solutions for WNM model for "Gal" 
dust. Curve labeled "P" is background generated by solid curve 
in Figure lb, in which optical pumping is included. Curve labeled 
"NP" is result withoug optical pumping. Figure 3c is the same as 
3b, except that the "SMC" star formation history of Figure 2b is 
used. 



respectively. The curves labeled "P" correspond to the 
star formation histories generated by the solid curves in 
Figures lb and 2b, where "P" indicates that optical pump- 
ing is included (where optical pumping is the mechanism 
by which the populations of the ground-term fine-structure 
states are mixed through UV excitations to higher lying 
levels [e.g. Sarazin et al. 1979]). The curves labeled "NP" 
show backgrounds generated by star-formation histories 
without optical pumping. As discussed in in Paper I the 
true WNM solution lies between these limits. The pre- 
dicted background for the "Gal" model is between 2.5 and 
5 times the 95 % confidence upper limit on Iebl, and 
between 3 and 10 times this limit for the "SMC" model. 
In both models, sources with z> Zmi„=1.6 generate back- 
ground intensities exceeding the 95 % confidence upper 
limit, indicating that DLAs in which C II* absorption 
arises in WNM gas produce more background radiation 
than observed. When the contribution of lower-z galaxies 
is included, the observed upper limits to the backgrounds 
are exceeded by much larger factors. Had we computed 
backgrounds generated by the dashed curves in Figures 
lb and 2b, we would reach the same conclusion; i.e., the 
WNM models are ruled out. 

This disproof of the WNM-dominated solution appears 
to be robust. The FUV mean intensities inferred for 4 
DLAs in which H2 is detected (Ge & Bechtold 1997; Sri- 
anand et al. 2000; Molaro et al. 2002; Levshakov et al. 
2002) are comparable to Go=1.7, the value found for the 
Galaxy ISM (Draine 1978; note Gq is Ji, in convenient 
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units of 10^^^ ergs cm"^ s~^ Hz^^ sr^^). Although these 
results may need to be corrected for suppression of opti- 
cal pumping (Sarazin et al. 1979), the implication is that 
<'!/'*> is significantly lower than required by the WNM 
models. Furthermore, the backgrounds predicted for the 
WNM models are conservative lower limits. This is be- 
cause we assumed the pressure of the two-phase medium, 
Peq , to exceed Pmin , the maximum pressure allowed for gas 
in a pure WNM phase. Values of Peg < Pmin would result 
in lower values of nwNM, hence higher p*. Note, that high 
values of are predicted even if C II* absorption arises in 
warm gas with temperatures below that predicted by our 
two-phase model. We considered scenarios in which C II* 
absorption arises in gas with T~1000 K; i.e., in thermally 
unstable gas like that predicted by Vazquez- Semadeni et 
al. (2000). Figure 3c in Paper I shows that at such tem- 
peratures the total cooling rate is considerably larger than 
the 158 /im emission rate, especially in the absence of opti- 
cal pumping. The resulting backgrounds are significantly 
above the 95 % confidence upper limit on Iebl when the 
emission from z < 1 galaxies is added to the contribu- 
tion from DLAs. Therefore, WNM models, or any model 
in which [C II] 158 /xm emission does not dominate the 
cooling rate, are unlikely to be correct. 

By contrast, CNM models in which < 'tp^,{z) > decreases 
with decreasing redshift are consistent with the background 
data. Consider models in which <V'*> equals a constant 
given by the average of all the inferred values of V* for 
the CNM models. These are depicted by dashed curves in 
Figure 3a and correspond to the "Gal" (G) and "SMC" 
(S) star formation histories shown as dashed curves in 
Figures la and 2a. Both models predict Iebl{z>0) to 
be above the 95 % confidence upper limit on bolomet- 
ric intensity. Moreover, Figures la and 2a show that 
at z=0 is significantly above the comoving SFR inferred 
from the luminosities of local galaxies. On the other hand, 
the solid curves in Figure la and 2a depict model star for- 
mation histories that are compatible with Pi.{z) inferred 
for galaxies at 2; < 1. These curves generate the back- 
grounds shown as solid curves in Figure 3a, which are 
compatible with the limits on Iebl- For these models 
< ijjt,(z) > at z = must be significantly below the « 
10-2.2 yj.-i ijp(.-2 determined at high redshifts both 
for "Gal" and "SMC" dust. At first, this appears to be 
inconsistent with the observation that deduced for 

local disk galaxies is comparable to 10^^ '' Mq yr^^ kpc"^ 
(Kennicutt 1998). But the latter star formation estimates 
were made by averaging over the corrected de'Vaucouleurs' 
radius, Rq, whereas DLA absorption at any redshift would 
occur out to an average radius, Rhi^'2Rq (Wolfe et al. 
1986). However, to compute p',, it is necessary to average 
tpt, over Rhi, not Rq. Therefore, the appropriate value of 
<■)/'*> for local disk galaxies should be reduced by a fac- 
tor of four to 10~'^ ° M0 yr~^ kpc~^. We conclude that if 
DLAs evolve into normal galaxies, their SFR per unit area 
has decreased significantly since z « 1.6. We emphasize 
this conclusion holds for averaged over R^i rather than 
Rq, which is normally used for computing in nearby 
galaxies (Kennicutt 1998) . ^ 

^When we compared the average Ic in DLAs to the average 



2.3. Consensus Star Formation Model 

By ruling out the WNM hypothesis, we have eliminated 
half the models discussed so far. Nevertheless, the remain- 
ing CNM models contain significant systematic uncertain- 
ties, as indicated by the scatter amongst p*{z) inferred 
from the various dust hypotheses (Figures la, 2a). Here 
we attempt to assess these errors as well as errors due to 
other effects, and to deduce consensus values for p*{z). 

To estimate the size of the systematic errors, we test 
the sensitivity of the CNM models to variations of cru- 
cial input parameters. We find p*{z) to be sensitive to 
changes in equilibrium pressure, Peq, and that the effect 
is similar in magnitude to the scatter in p'*(z) due to un- 
certainties in the composition and depletion level of dust. 
The results in Figures 1 and 2 were computed assuming 
Peq = {PminPrnaxY^'^ ■ Bccausc of thc Uncertainties sur- 
rounding this criterion (see discussion in Paper I), we now 
consider the effects of letting P^q vary between Pmin and 
Pmax- We find that decreases with increasing Pgq. As 
Peq rises, ucnm increases, which results in lower values 
of ijj^, for a fixed Ic (see Figure 5 in Paper I). Therefore, 
p* is a minimum when Peq=Pmax-i and a maximum when 
Peq =Pmin- We assumc the variances in are determined 
by differences between these limiting values of p, and thc 
means defined by Peq={PminPmaxY^'^ ■ It is possible to 
increase the variances in p* further by relaxing the stan- 
dard ratio of cosmic-ray ionization rate, CcR.i to the SFR 
per unit area, "0*: given in equation (9) in Paper I. But 
the consequent increase in p, is constrained by the upper 
limit on Iebl to logiop* < —0.2 M© yr^^ Mpc^"*, while 
the decrease is limited to logioP* > — 1-5 M© yr^^ Mpc~^ 
by thc observed ratios of C II to C I column densities, 
which become too small when (cR/tp* is more than twice 
the standard value (see § 5.1). Though it is possible for p* 
to attain these extreme values, it is more likely to remain 
within the standard limits, which justifies their use in con- 
structing the variances. Similar procedures were used to 
compute variances in p* due to variations in other input 
parameters such as the ratio of the radius to height, R/h, 
of the model uniform disk (sec Figure 4 in Paper I) . 

In order to compute consensus values of p*{z) in a given 
redshift bin, we considered all four possible combinations 
of "Gal" versus "SMC" dust composition and minimal ver- 
sus maximal depletion. For each of these four models we 
determined thc probability distribution of p'*(z) using the 
best-fit value and an error budget that included (a) vary- 
ing Peg from Pmin to Pmax, (b) varying the aspect ratio 
R/h between minimum and maximum values, and (c) the 
random errors appearing in Figures 1 and 2. Although 
we suspect that "SMC" dust and minimal depletion are 
more likely to be correct, we conservatively assumed that 
all four models were equally likely. We performed a Monte 

Ic for thc ISM (Paper I), wc found thc ratio, (Ic)dla/ (lc)iSM 
roughly equaled thc average dust-to-gas ratio of DLAs relative to 
that of the ISM, i.e. k. This indicated Go in DLAs to be roughly 
equal to the ISM value, G()=1.7, which corresponds to i/),=10— 
Mq yr^-' kpc~^. Does this contradict our finding that -(/>, at z=Q 
equals lO^"*'" Mq yr^-'^ kpc^'^? The answer is no bccausc thc high 
opacity of dust in the ISM indicates Go arises mainly from local 
sources within the optical radius, whereas thc lower opacity of dust 
in DLAs indicates Go inferred from C II* absorption lines is a global 
average, which includes the very low SFRs occurring outside the 
optical radius 



6 



STAR FORMATION IN DAMPED Lja SYSTEMS 



Carlo simulation drawing an equal number of simulated 



OfT 



1 -0.6 - 



n 



e 



o 

o 




Fig. 4. — Green data points depict p, and 68 % confidence errors 
for "consensus" model described in § 6.3. Dotted data points are 
galaxy data described in previous figures. Smooth curve is eyeball 
fit to "consensus" model at high z and galaxy data at low z. 

data points from each of the four model probability dis- 
tributions. Note that this is equivalent to a Monte Carlo 
simulation where each DLA is analyzed using all four mod- 
els and then these 4xn£)i^ data points are resampled at 
random to generate the maximum possible variance. The 
resulting probability distribution for /ci*(2;) is well described 
by a Gaussian, and we computed the resulting mean, 68% 
confidence intervals, and 95% confidence intervals. There 
are additional systematic uncertainties that we are unable 
to quantify at present, including those due to uncertainties 
in the grain size distribution, and others that we are unable 
to compute such as possible radial variations in -0* . We do 
not expect these additional sources of error to dominate. 
In particular, we show in § 3.2 that the error due to our 
assumption of spatially uniform tp^ is probably less than 
0.1 dex. The largest systematic uncertainty at present is 
produced primarily by the model with "Gal" dust compo- 
sition and maximal depletion, so falsifying cither the "Gal" 
or maximal depletion solutions would raise the result for 
p»{z) and significantly reduce the uncertainty. 

The consensus "Madau" diagram for DLAs is shown in 
Figure 4. The results are for the Einstein deSitter cos- 
mology, and the interested reader is referred to Table 1 
where we compare these with inferred for the standard 
A cosmology; the differences are of order 0.1 dex. The er- 
ror bars (corresponding to 68 % confidence levels) are of 
course larger than in Figures 1 and 2, which include only 
random errors. Though our errors are larger than reported 
for the Lyman-Break galaxies, the latter errors do not in- 
clude systematic errors such as extinction corrections to 
galaxy luminosities, which are surely present (see discus- 
sion in Steidel et al. 1999). By contrast the effects of dust 
are essential features of our models. Moreover the low val- 
ues of K imply that at least half of the radiation from the 
disk is emitted at rest-frame FUV wavelengths; i.e., our 



dust correction is less than a factor of 2. The blue curve 
is our eyeball fit through the DLA and low-^ galaxy data 
in this diagram and will be used in the following section 
to compute integrated quantities such as the mass in stars 
and metals produced over various time scales. Of course, 
it is possible that high-z DLAs do not evolve into low- 
z DLAs and their associated galaxies, but instead evolve 
into a population of objects with luminosity density far 
below that of normal galaxies. While we cannot rule out 
this scenario altogether, we believe it is implausible. The 
principal argument against it is the agreement between 
the comoving mass density of neutral gas in high-z DLAs 
and the mass density of visible matter in current galaxies 
(Storrie-Lombardi & Wolfe 2000). This indicates a con- 
nection between DLAs at high redshift with those at low 
redshift, unless one assumes this agreement is a random 
coincidence. As a result, the most likely scenario is one 
in which < 4'*{z) > decreases in time at z < 1.6 in such a 
way that the DLAs evolve into low-z galaxies. 

Figure 4 also shows consistency with p*(z) inferred from 
gas consumption in DLAs (Pei et al. 1999), which lends 
credibility to the idea that the decline with time of the 
comoving density of neutral gas is related to star forma- 
tion. Furthermore, Figure 4 indicates approximate agree- 
ment between p*(z) determined for DLAs and LBGs. That 
measurements of the same quantity by independent techniques 
based on different physical considerations are even in approximate 
agreement is cither a coincidence or indicates a connection 
between DLAs and LBGs. We shall address this issue in 

3. IMPLICATIONS 

In this section we discuss several consequences of this 
work. In particular we discuss (a) the production of stars 
and metals implied by the derived p'»(z), (b) a scenario 
in which star formation is confined to a centrally located 
bulge, and (c) evidence for feedback. 

3.1. Baryon and Metal Production 

Having determined /9*(z) for DLAs, we can integrate 
under the smooth curve in Figure 4 to obtain /0*(z), the 
comoving mass density of stars at redshift z. We find 

p*{z) — / p^,\dt/dz\dz , (7) 

J z 

where this expression for /3*(z) is independent of cosmol- 
ogy and Hubble constant when /0*(z) is determined from 

logioP* 

Moyr^^Mpc"-^ 
z Einstein dcSitter'^ Standard A" 
~^J5 -0.58±0.42 -0.68±0.42 
3.70 -0.59±0.35 -0.71±0.35 

Table 1 
Cosmology Dependence of p, 

"Cosmology with f2^=l, f2A=0, h=0.5 
''Cosmology with 0^=0.3, nA=0.7, h=0.7 
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the C II* technique; i.e., from equation (4). Note that un- 
hke Iebl, this integral receives considerable weight from 



was computed assuming Pmetais =(l/42)p* (Madau et al. 
1996), and is compared with the comoving density of 
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Fig. 5. — Mass density of stars per unit comoving volume for 
various populations. Solid blue horizontal line depicts density pro- 
duced by SFR for DLAs. Dotted blue lines are 68 % confidence level 
error bars. Red horizontal lines depict corresponding densities for 
stars formed by LBGs, current spheroids, current spiral bulges, and 
current spiral disks. "Gas in DLAs" corresponds to comoving mass 
density of gas in DLAs at 2; 3. 

z > 1.6. Pettini (1999) used the last equation to deter- 
mine p*(0), the current mass density in stars formed by 
the LBGs. Our estimate of p*(0) for DLAs, which we ob- 
tain by integrating the last equation from z—0 to zp, is 
shown in Figure 5 along with 68 % confidence intervals. 
We compare this to determinations of the current mass 
densities of stars formed by LBGs and other stellar popu- 
lations (Fukugita et al. 1998). In deriving this result we 
(a) integrated to the present under the solid curve shown 
in the figure, and (b) multiplied by Leitherer's (1998) cor- 
rection factor of 0.4 (adopted to correct for a more realistic 
IMF [see Pettini 1999]). Figure 5 shows DLAs and LBGs 
produce the same mass in stars to within 1 a. Moreover, 
the star formation history of DLAs suffices to produce 
the observed stellar content of spheroids, bulges of spi- 
rals, and spiral disks. Though the indicated uncertainties 
in p*(0) are large, the similarity between the predicted 
stellar content of DLAs and observed stellar content of 
galaxies is consistent with the idea that their progenitors 
were DLAs (e.g. Wolfe 1995). In addition, the similar- 
ity between the comoving gas density in DLAs at z ^ 
3 (Storrie-Lombardi & Wolfe 2000) and p*(0) is further 
evidence of self-consistency, though some infall might be 
required if the gas density is really lower than (0) . 

We also updated Pettini's (1999) calculation for the 
mass of metals produced by 2: = 2.5, the median redshift of 
DLAs for which metal abundances have been determined 
(Pettini et al. 1994; PW02). Pettini (1999) obtained this 
result using comoving SFRs for LBGs whereas we use co- 
moving SFRs for DLAs. The result shown in Figure 6 



Fig. 6. — Comoving mass density of metals. Solid blue line shows 
metals produced in DLAs by z = 2.5. Dotted blue lines are corre- 
sponding 68 % confidence contours, and dashed blue line is lower 
95 % contour. Red lines correspond to density of metals in spiral 
bulges and in DLAs. 

metals in 2; = 2.5 DLAs and with the current mass density 
of metals in spiral bulges. Clearly the metals produced are 
sufficient to account for the metal content of spiral bulges. 
However, as discussed by Pettini (1999) the mass density 
of produced metals is 30 times higher than metals observed 
in DLAs. The difference is significant, since as shown in 
Figure 6, the observed metal content of DLAs is well below 
the 95% confidence contour predicted for the produced 
metal content. Therefore, the difference between observed 
and produced metal content is real and leads to a "missing 
metals" problem. Pettini (1999) first noticed this problem 
when he found that the metals produced in LBGs exceeded 
the metals measured in DLAs. The problem is much more 
severe in our case because we are measuring both metal 
production and metal content in the same population. 

Three possible solutions to the missing metals problem 
are to sequester the metals produced away from the DLA 
gas observed at z 2.5, either by confining these met- 
als in "bulges" , in different systems, or in the IGM. The 
first solution requires that most of the star formation we 
are seeing is occurring in a compact region, i.e., a bulge, 
and metals produced in this region do not rapidly enrich 
the rest of the gas of the galaxy beyond the low metal- 
licities observed in DLAs. The second solution requires 
that DLAs are a transitory phase early in the formation of 
galaxies, meaning that by the time significant metals have 
been produced the neutral gas has already been used up. 
Hence the objects observed as DLAs ai z ^ 2.b are entirely 
distinct from objects that follow the DLA star formation 
history starting at formation redshifts z > 4. The third 
solution solution is to allow the metal-enriched material of 
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supernovae to blow out of the DLA galaxy (Mac Low & 
Ferrara 1999). We see two problems with blow out. First, 
the efficiency of ejecting metal-enriched gas must exceed 
1— (29/30)=0.97, which is much larger than the maximum 
efficiency of 0.5 seen in local starburst galaxies (Martin 
2003). Second, such ejection would result in a mean IGM 
metallicity, [M/H]=— 1.2, which is at least two orders of 
magnitude larger than the metallicity of the Lya forest 
(Songaila 2001); this would require placing most of the 
metals in the high— z IGM in some imdctcctable state. In 
the next subsection we shall explore the solution we find 
most appealing, the bulge hypothesis. Other possible ex- 
planations include changing the IMF from that assumed by 
Pettini (1999) (which is the same Madau et al. [1996] IMF 
used in our previous computations) so that lower masses 
of metals are released at the endpoints of stellar evolution. 

3.2. Bulge Hypothesis 

Suppose star formation is concentrated in the central re- 
gions of DLAs, e.g., in proto- bulges, which at early times 
could be configurations of molecular clouds that are rarely 
detected in DLAs because of a small covering factor or pos- 
sible obscuration by dust. In that case most of the metals 
would be released in the bulge, which would explain why 
the mass of metals produced in DLAs is consistent with 
the mass of metals in spiral bulges within the errors, as 
shown in Figure 6. The Milky Way bulge is relevant since 
it is a metal-rich but old population of stars most of which 
formed by z = 2.5 (Wyse et al. 1997). In this picture 
a small fraction of the metals produced would find their 
way to the outer disk via stellar winds or supernovae ex- 
plosions, thereby explaining the lower metallicities of the 
DLA gas. This idea is self-consistent since, as we now 
show, the higher star formation rate per imit area in the 
central region can account for the heating rates inferred 
from the C II* observations of the outer disk without sig- 
nificantly increasing p^,. 

For simplicity, let the bulge be a sphere with radius Rb, 
which is located at the center of a uniform disk with ra- 
dius R and half-thickness h. In this scenario the disk gives 
rise to damped Lya and C II* absorption, while FUV 
radiation emitted by the bulge is the source of the mean 
intensity at radius r, J^{r), which heats the gas through 
grain photoelectric emission. In that case the mean inten- 
sity is given by 

J?ir) = j^,eM-KT) , (8) 

where Rb < r < R, is the Luminosity per unit fre- 
quency bandwidth of FUV starlight emitted by the bulge, 
and ki, is the absorption opacity due to dust in the uniform 
disk. 

The radiation intensity inferred from DLAs by our method 
represents an average over all possible lines of sight through 
these uniform disks, where the C II* column determined 
for each system represents the average of all gas along that 
particular line of sight. We approximate this "average of 
averages" as a simple average of the radiation intensity 
received at all points in the disk, i.e. 
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with C = 2/[7ri?22M7r], g{r) = exp{-k,/r)/{4Trr'^), and 
where we have assumed Rb < h, which seems reason- 
able for the thick disks expected at high redshift. Ex- 
cept for the Rb limit, this expression is criuivalcnt to the 
solution for obtained for uniform disks (see eq. 14 
in Paper I) since the luminosity density of such disks, 
pi,=Li,/{'KR^2h). This occurs because of the symmetry 
between the flux received at the center of a uniform disk 
from all points within the disk and the average fiux re- 
ceived over a uniform disk from a bulge located at its cen- 
ter. When integrated, the solution looks just like that of 
equation 15 in Paper I except that the 1 is replaced by 
exp(— fc,yi?B), which will be very close to 1 given that the 
entire disk is nearly optically thin (the bulge of course is 
likely to be optically thick and have a different value of 
kv, but our observations are sensitive to the FUV pho- 
tons that successfully escape from the bulge region so this 
does not affect the results). If L^/{'KR^2h) is the same for 
the bulge and disk models, the expected values of radia- 
tion intensity are equal to an accuracy within 10 %, i.e. 

< Jjf >~ where is the mean intensity computed 
for disks with a uniform distribution of sources in equation 
(14) in Paper I. Since these two models represent extremes 
of the source distribution (uniform versus central source) 
we expect any intermediate source distributions to lead 
to similar values of radiation intensity. This is important 
since hierarchical structure formation implies that a given 
DLA could receive radiation from several compact regions 
of active star formation rather than a single central bulge. 
Because the SFRs per unit comoving volume for bulges or 
disks are given by 

cx LyUco , (10) 
the bulge-to-disk ratio {p*)b/{p*)d= L^/L^- Compari- 
son between equations 14 in Paper I and equation (9) in 
this paper shows this ratio equals the product of 

< >/ and the ratio of the dimensionless integrals 
in both equations. To compare bulge and disk SFRs we 
assume that J,^=< > so that both models generate 
the observed heating rate. Because the dimensionless in- 
tegrals are identical to within 10%, we find that {p*)b^ 
{p*)d- 

As a result the estimates of deduced for star forma- 
tion throughout uniform disks do not change significantly 
when star formation is confined to the centers of such disks. 
Therefore, our estimates of for disks do not appear to 
be very sensitive to the radial distribution of the sources 
of FUV radiation provided the disks are optically thin to 
such radiation. Though the disks are likely to be opti- 
cally thin, we cannot rule out the presence of optically 
thick dust in the bulge, which attenuates some fraction 
of the FUV radiation emitted by the stars. In that case 
the expression for (p,)b is found by equating p^, for the 
disk and bulge models and the resulting {p*)b would be a 
lower limit to the actual SFR per unit comoving volume. 
Since we can increase p*(^>1.6) inferred for the consensus 
model, i.e., inferred from FUV heating, by a factor of 3.3 
before the background is violated, the radiation attenu- 
ated in the bulge model can be as much as a factor of 2.3 
times that observed; i.e., as much as 0.7 of the total FUV 
radiation can be attenuated. 
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We can also use the bulge model to compute the total 
SFRs, M*, required to explain the observed C II* heating 
rates. From equation (8) we find that AI^— 1.9(r/10kpc)^Go 
M0yr~^, where we used the Madau et al. (1996) calibra- 
tion to convert L^, to M*. Assuming Go = 6.8, the average 
of the positive detections for the "Gal" minimal depletion 
model, we have M*=13(r/10kpc)'^ MQyr~^. If we let r 
equal 10 kpc, a typical impact parameter, we find that 
is consistent with upper limits obtained from Ha imaging 
of 7 DLAs by Bunker et al. (2001). On the other hand 
this SFR is about 3 times higher than the more sensitive 
NICMOS upper limits reported by Kulkarni et al. (2001) 
for the DLA toward Q1244+34. Unfortunately we do not 
have C II* 1335.7 profiles of this DLA. Future searches for 
emission from DLAs with detected C II* absorption may 
prove to be sensitive tests of the bulge hypothesis. These 
tests will be less sensitive for the uniform disk models, 
where the predicted surface brightnesses are low. In any 
case the limits on for the bulge hypothesis place severe 
constraints on rotating disk models in which r>20 kpc is 
required (PW97), since for the central bulge will be 
higher than observational constraints allow. 

3.3. Search for Evidence of Star Formation and Feedback 

Stars leave imprints on the gas from which they form, 
and this may be detectable in DLAs. In the uniform disk 
scenario a correlation should exist between -0* and 
iV(II I) resembling the Kennicutt relation found in nearby 
disk galaxies (Kennicutt 1998). This is a manifestation of 
the condensation of gas into stars, since it is equivalent to 
the Jeans instability criterion in rotating disks (Toomre 
1964). Stars also deposit energy and other byproducts of 
stellar evolution into the gas, and so one might expect to 
find evidence of feedback if the uniform disk scenario is 
correct. This would occur through shocks generated by 
supernova explosions. Because of the short main sequence 
lifetimes of the progenitor stars, feedback is directly re- 
lated to the SFR. By contrast, in the bulge scenario the 
connection between stars and DLA gas may (a) not exist 
or (b) be indirect. 

We test these ideas with SFRs derived from the CNM 
model with "Gal" dust and minimal depletion. The results 
of this investigation are statistically indistinguishable from 
those of "Gal" dust with maximal depletion, and "SMC" 
dust with maximal or minimal depletion. The test of the 
Kennicutt relation is shown in Figure 7a where we plot 
(iV(II I), -0*) pairs for DLAs along with the l-cr contours 
of the Kennicutt relation.^ Clearly the DLAs show no ev- 
idence for a Kennicutt relation. Specifically, a Kendall t 
test using only positive detections reveals t — 0.25 and 
the probability that the null hypothesis of no correlation 
is correct is p Kendall— 0.28. This conclusion differs from 
the findings of Lanzetta et al. (2002) who used the Ken- 
nicutt relation to deduce the SFR intensity distribution 
as a function of SFR per unit area from the frequency 
distribution of H I column densities in DLAs. They found 
excellent agreement with the intensity distribution deduced 
directly from galaxy brightnesses. In that case is de- 

Although the total H column density is used for nearby galaxies, 
A'^(H I) is adequate for DLAs where the molecular content of the gas 
is low (Petitjean et al. 2002); Ledoux et al. 2002 



1 1 1 i i 1 1 1 
(a) 


1 1 1 1 1 1 1 1 1 


(b) • . 


T 






|-L|'t '~t 


^ : 














20 20.5 21 21.5 -3 -2 -1 



log,oN(HI) [cm-s] [Si/H] 

Fig. 7. — (a) Comparison between A'^(H I) and ip* 
for DLAs with Kennicutt (1998) relationship defined by 
i/i. = (2.5±0.7xlO-'')x(Ar(H I)/l. 26x10^0 cm-2)i-4±0-i5 
Mq yr~^ kpc~-^. Data points inferred with "Gal" dust and 
minimal depletion. Red data points are positive detections, green 
are 2-cr upper limits, and blue is 2-cr lower limit, (b) Comparison 
between [Si/H] and i/i*. Same color coding as 7a. Magenta circle is 
estimated mean for LBGs. 



duced from emission from pixels with linear dimension of ~ 
0.2h~^ kpc. By contrast, the H I column densities in DLAs 
are sampled over transverse distances determined by the 
linear scale of the continuum source in QSOs, which is typ- 
ically less than 1 pc. As a result, the agreement between 
column density measurements on small scales with SFRs 
per unit area on larger scales indicates that on average 
the Kennicutt relation for nearby galaxies holds in high-z 
DLAs (Lanzetta et al. 2002). 

The reasons why the points in Figure 7a display so 
much scatter about the Kennicutt relation are straight- 
forward. First, if star formation occurs in the DLA gas, 
tp^ inferred from the C II* technique is averaged over the 
linear dimensions of the DLA, which exceed 5 kpc in any 
model (e.g. Haehnelt et al. 1998). Because N{li I) likely 
varies on scales smaller than 5 kpc, correlations between 
(7V(II I),?/;*) pairs are not expected. Second, in the bulge 
model, star formation does not occur in the gas giving 
rise to damped Lya absorption. As a result no correla- 
tions between (iV(II l),tp*) pairs are predicted. However, 
it may be possible to distinguish between the two mod- 
els with a sufficiently large data set. If star formation 
occurs in the DLA gas, 7V(H I) averaged over the DLA 
sample should correspond to iV(H I) averaged across a 
typical DLA, and therefore the sample averages of 7V(H I) 
and ■0* should obey the Kennicutt relation. On the other 
hand, no such correlation is predicted for the bulge model. 
Interestingly, the averages over the positive detections in 
Figure 7a result in <0*>=6.4x 10""^ Mq yr~^ kpc~^ and 
< A^(HI) >=1.0xlO^^ cm~^, which are within la of the 
average Kennicutt relation. In the bulge model this is 
merely a coincidence. 
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Fig. 8. — (a) Comparison between C IV velocity width Avciv 
and . (b) Comparison between low-ion velocity width Avi^^w a-nd 
tpt, . Color coding same as in previous figures. 



To search for evidence of feedback, we first looked for 
correlations between SFR per unit area and metallicity. 
Nearby spirals exhibit negative radial gradients in metal- 
licity (Garnett et al. 1997) and in SFR per unit area (Do- 
pita & Ryder 1994), implying a correlation between metal- 
licity and SFR per unit area. Such metallicity gradients 
may have also been detected in DLAs (Wolfe & Prochaska 
1998). We used ([Si/H], ip^) pairs to search for such cor- 
relations. The results shown in Figure 7b yield tentative 
evidence for a correlation since t=0.46 and PKendaU=Q-05, 
where again we used only positive detections. We then fo- 
cused on kinematic evidence for feedback. In the ISM, 
enhancements in velocity width are found in regions of 
higher than average SFRs such as Orion (Cowie et al. 
1979) and Carina (Walborn et al. 1998; Savage et al. 
2001). Evidently, gas in these regions is stirred up by in- 
creased supernova activity. We tested for kinematic feed- 
back by looking for correlations between absorption-line 
velocity width and V*- Specifically, we checked for kine- 
matic feedback in the neutral gas by searching for corre- 
lations between {Avlow, ip*) pairs, where Avlow is the 
velocity width of low ions in DLAs (Prochaska & Wolfe 
1997; hereafter PW97). We also checked for kinematic 
feedback in the ionized gas by searching for correlations 
between (AvciVi'ip*) pairs, where Avciv is the velocity 
width of the C IV 1550 transition. The results using only 
positive detections are shown in Figure 8a and reveal no 
evidence for correlations in the ionized gas, since t=0.97 
and PKendaii—0-76- But Figure 8b does reveal possible ev- 
idence for kinematic correlations in the neutral gas, since 
r=0.77 and p Kendall =0-002. 

The reasons for null correlations in the case of feedback 
are the same as discussed above. Namely, no correlations 
are expected for the bulge hypothesis because star for- 
mation occurs in regions disconnected from the absorbing 
gas used to infer metallicity and kinematics. In the disk 



model, ■0* is averaged over linear scales large compared 
to the ~ 1 pc transverse dimensions sampled in absorp- 
tion. Therefore, in this model the null correlation between 
Avciv and ip^, indicates significant random variations in 
high-ion velocities on linear scales small compared to the 
length scale of the star forming regions in DLAs. Because 
it is reasonable to expect similar variations in the case of 
low-ion velocities, how can we understand a correlation 
between the {Avlow, V'*) pairs, if confirmed? The answer 
may be that Avlow has global rather than local signifi- 
cance. Specifically, Avlow may reflect the depth of the 
gravitational potential well of the DLA, as predicted in 
the case of rotating disks (PW97) or protogalactic clumps 
(Haehnelt et al. 1998); i.e., the SFR per unit area may be 
correlated with total mass. Note, this explanation would 
apply both to the uniform disk and bulge models. A simi- 
lar explanation might also apply if the tentative correlation 
between ([Si/H],';/)*) pairs is confirmed; namely, that the 
metallicity of the gas does not vary randomly on scales 
small compared to scale of the galaxy hosting the DLA. 
Rather, metallicity might be a function of total mass, as 
in the case of current galaxies. 

Of course all these results need to be tested with more 
data. In particular, the statistical significance of the cor- 
relation between Avlow and -0* would be reduced if the 
two upper limits with > 150 km s^^ were added to the 
sample of positive detections. On the other hand, the two 
DLAs with logio0*«— 3.6 yr~^ kpc~^ are the "out- 
liers" discussed in Paper I. These are more likely to be 
WNM-dominated absorbers with significantly higher val- 
ues of 0*, in agreement with the predicted correlations. 
4. TESTS OF THE MODELS 

Having presented evidence for star formation in DLAs, 
and having the described the implications of the derived 
star formation histories, we now discuss three tests of the 
models upon which these results are based. The first is 
a test for the grain photoelectric heating mechanism, the 
second tests the hypothesis that Ic is a cooling rate, and 
the third describes a search for evidence of a two-phase 
medium. 

4.1. Evidence For Grain Photoeleetrie Heating 

A critical premise of the C II* technique is that neutral 
gas in DLAs is heated by photoelectrons ejected from inter- 
stellar grains by FUV radiation emitted by massive stars. 
In that case the heating rate per H atom. Yd oc Ke'i/'*- The 
efficiency of grain photoelectric heating, e, is essentially 
constant as it is insensitive to variations in electron den- 
sity, temperature, and FUV radiation field in the CNM. 
Moreover, the scatter in -0* is limited such that the av- 
erage logio<'0*>=-2.19lS];26 Mq yr-i kpc~^ Therefore, 
a prediction of the grain photoelectric heating scenario is 
that Yd should be roughly correlated with the dust-to-gas 
ratio, K. 

We test the grain photoelectric heating hypothesis by 
letting lc=Yd (which is an excellent approximation in the 
CNM) and then compare 1^ with k. The two quantites are 
plotted against each other in Figure 9. The red data points 
are positive detections, green are upper limits, and blue is 
a lower limit. The solid curves are lines of constant pre- 
dicted by a version of the CNM model with a fixed redshift 
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Fig. 9. — Plot of dust-to-gas ratio, re versus Ic where the "Gal" 
dust and minimal depletion is assumed to compute re. Color code 
of data points same as in previous figures. Continuous curves are re 
versus Icr relations for constant logio^/)* = —4.0, —3.0, —2.0, —1.0, 
and 0.0 Mq yr~^ kpc~^ predicted by model described in text. 



and a given prescription for computing n from metallic- 
ity. That is, for a given metallicity and '0* we calculate 
two-phase equilibria of gas subjected to grain photoelec- 
tric heating and assume the DLA density to be given by 
the computed ucnm- We include the CMB contribution 
to radiative excitations of the [C II] fine-structure states 
by assuming z — 2.8, the median redshift of the sample, 
and ignore optical pumping (which should lead to no loss in 
generality since optical pumping is negligible in the CNM) . 
To calculate k we adopt the CNM "Gal" model with min- 
imal depletion. As a result we let the intrinsic carbon 
abundance, [C/II]j„t = [Si/II]i„t-|-[Fe/Si]i„t and [Fe/Si]i„f 
= —0.2. We also assume [Fe/Si]gas=— 0.4, which is the 
average value for our sample. We then compute Icr with 
techniques described in Paper I. Visual inspection of Fig- 
ure 9 indicates several phenomena. First, the measured Ic 
and K are correlated. Specifically, performing the Kendall 
tau test using the positive detections alone we find t—O.QA 
and PKendaii^O -0064:, where pxendaii is the probability of 
the null hypothesis of no correlation. Second, the slope of 
the data in the (k, ^c) plane is approximately parallel to 
the model predictions, and all the positive detections are 
bounded by logioV'* = —2 and —3 yr^^ kpc"^. Third, 
the case for correlation receives additional support from 
the location of the upper limits at low dust-to-gas ratios, 
i.e., logioK < —2.2, and the lower limit at the relatively 
high dust-to-gas ratio, logioK > —1.6. 

Can the correlation between Ic and k be explained by 
heating mechanisms other than grain photoelectric heat- 
ing? The obvious alternatives are cosmic-ray and X-ray 
heating. Suppose the actual SFRs are orders of magni- 



tude lower than we infer for grain photoelectric heating, 
but the cosmic-ray ionization rate, QcRi is large; i.e., the 
ratio CcbJiP* is orders of magnitude larger than given in 
equation (9) in Paper I. In that case it is possible for the 
cosmic-ray heating rate, Tcr, to dominate the heating rate 
in the CNM. During ionization by cosmic rays, primary 
and secondary electrons are mainly liberated from H and 
He (W95). As a result, Tcr will be independent of metal 
abundance, hence independent of k. Since we have linked 
Cc_R, hence Tcr, to the SFR, and Ic~Tcr when cosmic 
rays dominate, Ic will be independent of k for a fixed -(/;*, 
in contrast to grain photo-electric heating. 

On the other hand, the X-ray heating rate could depend 
on metallicity, hence on k, since abundant heavy elements 
dominate the X-ray photoionization cross-section per H 
atom at photon energies above the Oxygen edge at 0.53 
keV (Morrison & McCammon 1983). We used the X-ray 
heating model of W95, which consists of two local thermal 
sources with T ~ 10^ K and an extragalactic power-law 
component. Heating by the thermal sources is dominated 
by photoionization of H and He since their X-ray spectra 
cutoff below 0.53 keV. Although the power-law component 
extends to energies above the Oxygen edge. X-ray heating 
in this case will also be dominated by photoionization of H 
and He. The reason is that for typical velocity-component 
column densities, iV(H I) « 1x10^° cm~^, most of the X- 
rays penetrating the H I gas have energies below 0.53 keV. 
Photoionization of H and He by these X-rays dominates 
photoionization of heavy elements by higher energy X-rays 
due to the low metallicities of most DLAs and the shape of 
the power-law spectrum. Thus, photoionization of heavy 
elements will not be the dominant source of primary elec- 
trons. Consequently, while the condition Ic ~ Txr might 
hold in the limit of low SFRs and high X-ray luminosities, 
Tx R. will be independent of k for a fixed tp* ■ 

To conclude, the correlation between Ic and k is natu- 
rally explained by grain photoelectric heating. The cor- 
relation at a fixed SFR follows from the physics of grain 
photoelectric heating, while the scatter of Ic at fixed k 
reflects the frequency distribution of V'*- Note, this corre- 
lation does not distinguish between the uniform disk and 
bulge models. Because the heating rate is the product of 
a "global" quantity, ip^,, and a "local" quantity, k, clouds 
with larger k will have higher heating rates independent of 
whether the incident FUV radiation arises locally or from 
the distant bulge. By contrast, neither the cosmic-ray nor 
X-ray heating rates are correlated with k at fixed SFR. 
To explain the observed correlation, one must postulate a 
correlation between -0* and k (or possibly metallicity). It 
is difficult to understand the origins of a correlation be- 
tween the "global" 0* and the "local" dust-to-gas ratio in 
this case. 

4.2. Is C ir Excitation Due to the CMB? 

Our models predict that in the CNM of most DLAs the 
cooling rate, nA, equals the spontaneous emission rate of 
158 /im radiation, Ic- This condition, nAfulc, holds when 
the ^P3/2 and fine-structure states in the ground 

term of C"*" are populated by coUisional rather than 
radiative excitations. In Paper I we described how CMB 
radiation populates these states directly, and how the FUV 
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Fig. 10. — Plot of Ic versus redshift. Red, green, and blue filled 
circles have usual meanings. Magenta stars depict {lc)cMB corre- 
sponding to measurements of Ic at each redshift. The plot shows 
data for 30 DLAs. Computation of {lc)cMB explained in text. 



radiation field, Gq, populates them indirectly through op- 
tical pumping via higher energy states. While optical 
pumping is important in the WNM, it can be neglected 
in the CNM. Although an increase in Go drives up the 
pumping rate, it also increases the coUisional excitation 
rate. This is because an increase in Go, and thus "0*: in- 
creases the grain photoelectric heating rate, which raises 
the CNM density, ucnm, as shown in Figure 5 in Paper I. 
As a result, the ratio of coUisional to optical pumping rates 
always exceeds unity in the CNM. ^ This need not be true 
for the ratio of coUisional to CMB excitation rates because 
the CMB does not heat the CNM gas (through Compton 
heating) in the redshift range of our sample DLAs. Be- 
cause the value of ucnm does not rise with increasing 
CMB intensity, the CMB may dominate collisions as a 
source of excitation in DLAs with low values of ucnm 
and high redshifts where the CMB intensity is high. 

Because the balance between coUisional and CMB ex- 
citations depends on various assumptions included in our 
models, it is important to assess their relative importance 
with model-independent tests. Figure 10 illustrates the 
results of a test relying on one free parameter, the carbon 
abundance, (C/H). Here we plot Ic, inferred from mea- 
surements of N{G II*), versus redshift for a sample of 
30 DLAs. For each DLA we also plot, as magenta stars, 
(IcjcMB versus redshift where {Ic)cmb is the value as- 
sumed by Icr when the CMB is the only source of excitation 
and de-excitation (in Paper I we showed that {IcrjcMB — 



2(C/II)A„;/li^„;exp 



CMB 



where the 



^Although the density of the WNM, nvi'jvM) a^lso increases with 
Go, the densities are not high enough for collisions to dominate op- 
tical pumping in most cases. 



present CMB temperature, Tcmb =2.728 K. Twenty five 
of the DLAs are drawn from the minimal depletion sample. 
Here we calculate {Ic)cmb by assuming [C/H] = [Si/H]— 0.2 
to compute (C/H). The plot also shows data for 5 addi- 
tional DLAs. They include the DLAs toward Q0201+11 
and Q2344+12 for which we assumed maximal depiction 
to obtain (C/H) because [Fe/Si] > -0.2 in both DLAs. 
We also included the DLAs toward Q0951-04, Q1425+60, 
and Q1443-f27 for which we assumed [C/H] = [Si/H]-0.2. 
These objects were excluded from the original minimal 
depletion sample because observational limits on [Si/H] or 
[Fe/H] prevented us from computing the dust-to-gas ratio, 
K, which is not needed to compute {Ic)cmb- The figure 
demonstrates that the CMB cannot explain the observed 
level of C II* excitation for the following reasons: First, de- 
spite the dispersion in carbon abundance, {Ic)cmb should 
increase rapidly with increasing redshift whereas the ob- 
served values of Ic show no dependence on redshift in the 
interval z = [1.6, 4.5]. Second, with two exceptions, the 
CMB excitation rate is too low to explain the observations. 
This is true even at z > 3 where the predicted values of 
{lc)cMB merge with the observed spontaneous emission 
rates. 

The exceptions are the DLAs at z=3.608 and 0=4.080 
toward Ql 108-07 and Q2237-06, the "outUers" discussed 
in Paper I. In these DLAs, {Ic)cmb comprises a signif- 
icantly higher fraction (> 0.8) of Ic than for the other 
DLAs. Because the Ic's are upper limits, the data place 
lower than usual upper limits on the [C II] 158 /im cool- 
ing rates. Wc interpreted this to mean that these sightlines 
pass only through WNM gas subjected to SFRs within the 
range determined from the CNM hypothesis. However, we 
cannot rule out the possibility that the gas is subjected to 
negligible SFRs, which lead to gas densities so low that col- 
lisional excitations are unimportant. In that event, CMB 
excitation alone would be responsible for the observed Ic- 
To decide between these hypotheses, we need to measure 
Ic and {lc)cMB more accurately. Though Ic can be deter- 
mined more precisely through better measurements of the 
C II* 1335.7 velocity profiles, higher signal-to-noise will 
not improve the accuracy of {Ic)cmb where the only source 
of error is in determining (C/H). The carbon abundance 
is difficult to obtain directly because the principal C II 
transitions are always saturated (e.g. Dessauges-Zavadsky 
et al. 2001). Instead we compute {Ic)cmb with the min- 
imal depletion assumption, [C/H] = [Si/H] -I- [Fe/Si] i„t and 
[Fe/Si]i„t=— 0.2, for 28 of the 30 data points shown in 
Figure 11. For the remaining 2 DLAs, [Fe/Si]gas > -0.2, 
which violates the nucleosynthetic ceiling limit in the min- 
imal depletion model, and we thus assume [Fe/Si]i„(=0; 
i.e., maximal depletion. We favor minimal depletion for 
most DLAs because it leads to self-consistent determina- 
tions of {lc)cMB, whereas maximal depletion leads to the 
inconsistent condition {Ic)cmb > Ic for the Q1108— 07 and 
Q2237— 06 DLAs. In any case the limited accuracy for de- 
termining (C/H) results in errors in {Ic)cmb of order 0.2 
dex, which is inadequate for distinguishing coUisional ex- 
citation from CMB excitation. 

In our treatment of CMB excitations of the C II fine- 
structure states we assume the CMB temperature at red- 
shift, z, Tcmb{z) — {1 + z)TcMB, whereas most published 
test this assumption by attempting to measure Tcmb{z) 
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Table 2 

PHYSICAL PARAMETERS FOR DLA GAS 





^1232^08 


goal 


7-38 


Parameter 


Estimate"^ 


CNM Model" 


Estimate'^ 


CNM Model" 


rt (cm ^•^) 


20-335 


7-14 


4-14 


3-6 


T(K) 


85-285 


130-188 


050'' 


130-200 


ne(cm 3) 


0.02 


0.002-0.03 




0.002-0.014 


Go 


~1.7= 


6.0-21.5 


~ 1.7= 


2.7-11.0 



"Prom Srianand et al. (2000). 

''Parameter range predicted by CNM models with "Gal" and "SMC" dust and minimal and maximal depletion. 
=Prom Levshakov et al. (2002). 

"^Uppcr limit from velocity dispersion a = 1.4 km s^^ observed for H2 lines differs from estimate of Levshakov et al. 2002. 
''Inferred from reported photoabsorption rate in Lyman and Werner bands, /3o 



directly. In principle this can be done by first making 
independent estimates of the collisional contribution to 
the level populations of the and ^P\/s fine-structure 

states, and then computing the black-body temperature 
required to explain the observed population ratio (e.g. Mo- 
laro et al. 2001). Accurate UVES echelle spectra were ac- 
quired by the VLT for two DLAs in our sample; Q0347— 38 
(Levshakov et al. 2001; Molaro et al. 2001) and Q1232-I-08 
(Srianand et al. 2000). In Table 2 we compare physical 
parameters obtained by these authors with predictions of 
our CNM models. The estimates are independent of our 
predictions because physical models of the DLA gas were 
not constructed. Rather techniques such as measuring 
velocity line widths to estimate the kinetic temperature 
or inferring Go from the fractional abundance of H2 were 
used. In some cases relative abundances of various rota- 
tional levels of H2 were used to infer the temperature, and 
the ratio Mg+/Mg° was combined with Go and photoion- 
ization equilibrium to obtain rig. The results in Table 2 
show, with the possible exception of Go , reasonable agree- 
ment between our predictions and these estimates. They 
are also inconsistent with the WNM hypothesis for both 
DLAs. 

4.3. Ratios of C IF to Resonance-line Velocity Profiles: 
Probing the Two-Phase Medium 

Here we discuss an observational test of the two-phase 
medium, a key clement in our models of DLAs. We de- 
scribe various aspects of the test and leave more quan- 
titative analyses for future discussions. Note, Lane et 
al. (2000) provide independent evidence for a two-phase 
medium from their analysis of 21 cm absorption in a DLA 
with z = 0.0912. 

Suppose, by analogy with the ISM, the probabilities 
that the line of sight intercepts WNM and CNM clouds 
are comparable (Kulkarni & Heiles 1987). If velocity com- 
ponents (i.e., clouds) in each phase have comparable H I 
column densities and the same clement abundances, then a 
random sightline through a typical DLA should encounter 
similar column densities of metals in each phase. This can 
lead to measurable differences between the velocity pro- 
files of C II* 1335.7 and resonance transitions such as Si 
II 1808.0, as we now show. 

Consider a two-phase configuration in which the frac- 
tions ycNM and ywNM of the total A''(H I) are in the 



CNM and WNM, where ycMN+ywNM=i- Assume the 
systemic cloud velocities to be ucnm and uwnm, and the 
internal Gaussian velocity dispersions to be ucnm and 
(^WNM- In that case the optical depths at velocity v in 
C II* 1335.7 and Si II 1808.0 are given by: 



ire 



1 



/cii* Acir 



{la 



yCNM 



$(u-MCiVM) 



mc J2n 



jV(HI) < 



I ilc)wNM ^(v — UwNm) 



/siiiAsiii(Si/H) 



yCNM 



CII* 



+yWNM 



•S>(v-uwnm) 



Sill 



where we used equation (3) in Paper I to compute 
r„(C II*), the /'s and A's are oscillator strengths and tran- 
sition wavelengths, {Ic)cnm and {Ic)wnm are Ic in the 
CNM and WNM, $(u) is the velocity profile normalized 
such that / ^{v)d{v / y/2T:a)=l, and (Si/H) is the Si abun- 
dance. 

To compute {Ic)cnm and {Ic)wnm, we evaluate Ic at 
the phase densities ncNM and uwnm- Equation (11) in 
Paper I shows that Ic = nAcii + {lc)pump + {Ic)cmb where 
nAcii is the net loss of thermal energy due to 158 /im 
emission and {lc)pump and {Ic)cmb are the spontaneous 
energy emission rates in the limits of pure optical pumping 
and CMB excitation, which are defined in equation (12) 
in Paper I. At the high values of ucnm, the fine-structure 
states are mainly populated by collisions and as a result 
Ic equals the cooling rate nh.cn, which equals the grain 
photoelectric heating rate in the CNM; i.e., 

{lc)cNM = Ld . (12) 

But at the low values of nwNM any one of the of the 
3 terms in the above equation for Ic can dominate. For 
example {Ic)wnm~{Ic)cmb in the 2;=3.736 DLA toward 
Q1346-03 for the CNM solution in which Go=1.95. By 
contrast {lc)wNM^{lc)pump in the z=2.039 DLA toward 
Q0458-02 for the CNM solution in which Go>13.5 (see 
Figure 5 in Paper I). Without optical pumping, {Ic)wnm 
~nAcii if the metallicity is low and observed Ic is high. 
In this case Ic does not equal because [C II] emission is 
not the dominant coolant. If {Ic)wnm is governed by op- 
tical pumping, the ratio {Ic)wnm / {Ic)cnm is rather well 
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determined because the quantities and T^i are both 
proportional to the radiation intensity Go and to metalHc- 
ity. Combining the last equation with equations (1) and 
(7) in Paper I and assuming [C/H]gas = [Si/H] + [Fe/Si]i„t 
we find 

ilc)wNM ^ „ c^irig jC/B^Q (Tiu) ISM 

{lc)cNM ~ " [l-10[Fc/Si],.e-[Fc/Si]„.],^^^,^ ' 

where {Tiu)isM, the optical pumping rate corresponding 
to the Draine (1978) radiation field, equals 1.57x 10"^° s"^ 
(Silva & Viegas 2002), and ccnm is the heating efficiency 
in the CNM. In deriving the last equation we assume the 
pumping rate corresponding to intensity Go is given by 
r;„(Go) = (r;„)/5M(Go/1.7). Averaging over the positive 
detections in the "Gal" minimal depletion model we find 
the average ratio <{Ic)wnm / {Ic)cnm> =0.097±0. 020 with 
optical pumping and <{Ic)wnm / {Ic)cnm> =0.067±0.011 
without pumping. As a result, equation (11) implies the 
WNM is undetectable in C IF 1335.7 for any value of 
ywNM, but is detectable in resonance transitions such as 
Si II 1808.0 for ywNM > 0.1. Consequently, a significant 
decrease in rt,(C II*)/T„(Si II) from the mean would signify 
the presence of the WNM. 

Such variations may have been detected in the DLA to- 
ward Q0347— 38. Though Figure 11a shows clear evidence 
for distinct velocity components at v——8 km s^^ and 
^;=+12 km s-i in Si II 1808 and Fe II 1608, the v=~8 
km s~^ component is not seen in C II* 1335.7. Rather, 
an asymmetric blue wing extends from v=+12 km s^^ to 
v=—20 kms. Because some of this absorption is due to 
weak C II* 1335.66 (/i335.66//i335.7i=0.11), the observed 
optical depth of the wing places an upper limit on Ty{C II* 
1335.71). Therefore, we conjecture that the v=—8 km s^^ 
component consists of WNM gas, while the i;=+12 km s^^ 
component consists of CNM gas. The presence of H2 ab- 
sorption only at t;=-|-12 km s"""^ and the lower S"'"/N° ratio 
in this component (Levshakov et al. 2002) supports this 
interpretation. Notice that the relative optical depths of 
these two components in Fe II 1608 is different than in 
Si II 1808. The most likely explanation is enhanced Fe 
depletion in the v=+12 km s~^ component. This does 
not affect our interpretation of the C II* profile because 
enhanced depletion at w=-|-12 km s~^ cannot explain the 
weakness of C II* 1335.7 at v=—8 km s^^ . In any case the 
significant decrease in the ratio, TZ{v)= Ty{C II*)/ri,(SiII), 
with decreasing velocity shown in panel 4 of Figure 11a is 
naturally explained by the presence of WNM gas at nega- 
tive velocities. 

Figure lib compares C II* and resonance- line profiles 
for the DLA toward QOlOO+13 (PHL 957). By contrast 
with Figure 11a, the C II* profile in this case closely re- 
sembles the Si II 1808 and Ni II 1741 profiles. The only 
significant difference is at v <— 8 km s~^ where C II* 
1335.71 exhibits a blue asymmetric wing that is missing 
from Si II 1808 and Ni II 1741. We believe this is a blend 
with weak v = absorption in C II* 1335.66. Because we 
ignore the increased C II* absorption at v >-|-30 km s^^ 
, which is likely to be a blend with Lya forest absorp- 
tion lines, we find no compelling evidence for a WNM in 
this DLA. However, this does not rule out the presence 
of a WNM. Suppose the velocity components at u « -1-3 
km s"^ and -|-15 km s^^ each contain WNM and CNM 
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Fig. 11. — Comparison between C II* and resonance-line veloc- 
ity profiles for 3 DLAs. Top row of panels show profiles for C II* 
1335.7. Second row shows profiles for Si II 1808, and third row 
for Fe II 1608 (a) and Ni II 1741 (b and c). Bottom row shows 
logio [(7^— <7^.>)/<7?.>] computed for successive 20 km bins 
where TZ is ratio of C II* 1335.7 to resonance-line optical depth. 
Columns (a), (b), and (c) show data for DLAs toward Q0347— 08, 
QOlOO+13, and Q2231-00. 



gas with the same velocity dispersion. From equation (11) 
we see the ratio 7?,(w) will not vary with v. But, in a sce- 
nario where the multi-phase structure resembles that of 
the ISM, the Si II velocity profiles would be wider than 
the C II* profile, since awNM^^acNM in the ISM (Kulka- 
rni & Heiles 1987). In principle this is a signature of the 
WNM. But differences between the observed profiles will 
be diluted by the finite resolution of HIRES (FWHM«7-8 
km s"'^ ), which is comparable to the widths of most com- 
ponents (i.e., they are unresolved), and Poisson noise. A 
detailed evaluation of these effects will be discussed else- 
where. 

Figure 11c compares C II* and resonance- line profiles 
for the DLA toward Q2231-00. Though the profiles are 
noisier in this case, some effects are clear. First, the ratio 
Tl{v) does not exhibit significant variations at w < —60 
km s~^ . On the other hand, there is evidence for varia- 
tions of TZ(v) in the component centered at ?; = 
km s^^ . At u « -I- 10 km s~^ , Tl{v) is lower than <TZ>, 
the mean TZ{v) integrated over the entire profile (see panel 
4 of Fig. 12c). This is consistent with WNM gas a.t v ^ + 
10 km s~^ . At the same time, there is tentative evidence 
for enhanced C II* absorption between v ^ — 30 km s~^ 
and — 10 km s^^ where C II* 1335.71 exhibits stronger 
absorption than either Si II 1808 or Ni II 1741. While weak 
C II* 1335.66 absorption at w w km s^^ may cause the 
excess C II* 1335.71 absorption at w > — 10 km s~^ it 
cannot explain the excess C II* 1335.71 absorption at w < 
— 10 km s^^ . 

Excess C II* absorption can be due to increased heating 
of the CNM. In the case of absorption by CNM gas alone 
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we have 



_ T.(cir) 

' r„(SiII) 



/sill^Sill hVulAul 



10- 



(Si/H) 

where we combined the definition of Yd (equation 1 in 
Paper I) and equation (12). Detectable variations in 7?,(w) 
in CNM gas are most hkely caused by variations in e and 
Go rather than in absolute element abundances. This is 
because 



|^20[Fe/Si], 



lOP-'/Si] 



(Si/H) 



(Si/H) 



(15) 



and the recent analysis by Prochaska (2002) showed evi- 
dence for remarkable uniformity in the relative abundances 
of DLAs (with the exception of the DLA toward Q0347-38). 
As a result, 



e(w)Go(w)- < eGo > 



(16) 



< eGo > 

where 5ti{v) = \R{v)-<TZ>]/ < TZ >, and e{v) and Go(w) 
are the grain photoelectric heating efficiency and FUV 
mean intensity at velocity v. Because e is a function of 
GoVr/ue (Bakes & Tielens 1994; Weingartner & Draine 
2001a), Sti(v) will vary with v if any one of these parame- 
ters changes with v. If GoVT/ue « 5x10'^ K^/^cm^, the 
grains are mainly neutral and e is insensitive to variation in 
Go. Therefore, in this CNM limit, 6n{v) « AGo/Gq. The 
excess C II* absorption in the Q2231— 00 DLA can then re- 
sult from the incidence of a larger-than-average radiation 
intensity, Go, on the CNM gas in the velocity range v — 
—30 to —10 km s^^ . Spatial variations in radiation inten- 
sity along the line of sight could result from (a) nearer than 
average displacement of OB stars from this gas, (b) differ- 
ing CNM cloud distances from a centrally-located bulge 
source, or (c) in the case of CDM galaxy formation sce- 
narios, passage of the line-of-sight through separate pro- 
togalactic clumps with independent SFRs (Haehnelt et al. 
1998). The CDM scenario should produce standard devia- 
tions, ~ 1, because the SFRs would be uncorrelated at 
the different clump locations. But in the bulge scenario, 
we find AGq/Gq ^ h/b where h is the thickness of the 
surrounding disk and b is the sightline impact parameter. 
This would lead to ct^^ < 1, for typical values of 6. To 
discriminate between these hypotheses, larger samples of 
DLAs with measured TZ{v) are required. 

5. DO H IGH-Z DLAS CONTAIN CNM GAS ? 

In this section we discuss possible objections to the pres- 
ence of CNM gas in DLAs. Recall that if C IF 1335.7 ab- 
sorption arises in WNM gas, the inferred FUV radiation 
contributes more background radiation than observed. On 
the other hand if it arises in CNM gas, the background ra- 
diation is consistent with observations. Nonetheless, the 
following arguments have been made against the presence 
of CNM gas in high-z DLAs. 

5.1. C II/C I Ratios 

The first of these has to do with the large values ob- 
served for the ratio, iV(C II)/iV(C I) (hereafter C II/C I). 
Liszt (2002) constructed two-phase models similar to ours 
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Fig. 12. — Plot of observed versus predicted C II/C I ratios for 19 
DLAs from Table 1. With the exception of the filled red circle, which 
is a positive detection, the observed C II/C I ratios are lower limits. 
The filled blue circles correspond to DLAs for which we have positive 
detections of Ic- Filled green circles correspond to DLAs with limits 
on Ic and filled magenta circle corresponds to DLA with lower limit 
on Ic (sec text), (a) Results for model with logio A''m=19.3 cm~^ 
and standard ratio of Ccr/''P* (equation 9 in Paper I), (b) Same 
as (a) except logioA'^u,=20 cm~^. (c) Same as (b) except Ccfl/V'* 
equals 0.3 times the standard ratio. 



and computed C II/C I ratios for 5 DLAs for which ob- 
served ratios were available. While the ratios he predicted 
for the CNM were lower than observed, those predicted for 
the WNM were consistent with observed values. He con- 
cluded the observed C II/C I ratios permit no more than 
a few percent of the DLA gas to be in the CNM phase. 

We repeated these calculations with the CNM model 
with "Gal" dust and minimal depletion. In Figure 12 we 
compare the results with C II/C I ratios deduced for 19 
DLAs in our sample. With the exception of the filled red 
circle depicting a positive detection, the observed C II/C I 
ratios are lower limits, while the predicted C II/C I ratios 
are definite numbers if Ic is detected (filled blue circles), 
and upper or lower limits if limits are placed on Ic (filled 
green circles; see caption to Figure 12). Figure 12a shows 
results for a model incorporating many of the assumptions 
made by Liszt (2002). In particular, we assume the inci- 
dent soft X-ray radiation is attenuated by gas with H I 
column density logioA^Mi = 19.3 cm~^. We are in general 
agreement with Liszt (2002) that this CNM model cannot 
account for the large values observed for C II/C I. Al- 
though the model predictions should lie to the right of the 
diagonal line representing equality between observed and 
predicted ratios, one third of the sample is on or to the left 
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of the line. This disagreement takes on significance when 
it is reahzed that all the discrepant points are lower limits 
on the observed C 11/C 1 ratio and hence the actual ra- 
tios are displaced even further from agreement with model 
predictions. 

However, the predicted C II/C I ratios are sensitive to 
changes in model input parameters. This is evident in Fig- 
ure f2b showing results for the same model as in Figure 
12a except N^, is increased to 10^° cm~^; i.e., to the at- 
tenuating column density, which is more appropriate for 
DLAs and which is used in all our calculations (see discus- 
sion in Paper I). Comparison with Figure 12a reveals an 
increase in predicted C II/C I ratios by ~ 0.4 dex with a 
consequent improvement between theory and observation. 
The reasons for the increase in C II/C I are as follows: 
The increase in TV^, results in a reduction in X-ray inten- 
sity, which in turn causes a decrease in the heating rate. 
This decreases the gas pressure, which brings about a de- 
crease in the CNM density, ucnm as illustrated in Figure 
5 in Paper I. For a given Ic a decrease in ucnm causes an 
increase in "0* (hence Go), a decrease in Tig, and a rise in 
T . All three factors conspire to increase the C II/C I ratio, 
since it is proportional to Go/(fT.eaci(r)), where aci{T), 
the case A recombination coefficient to C I, decreases with 
increasing T. For many DLAs even better agreement is 
achieved if we retain logioA^u, = 20 cm^^ and reduce the 
ratio of cosmic ray ionization rate to SFR per unit area, 
CcrI^*, below the ratio given in equation (9) in Paper I. 
At these low X-ray intensities, cosmic rays still dominate 
the ionization rate, and as a result tlcnm is reduced even 
further. Figure 12c shows results for CcR/i^* equaling 0.3 
times the ratio in equation (9) in Paper I. In this case most 
of the predicted C II/C I ratios are in better agreement 
with observations than before. However, in the most metal 
rich DLAs the additional decrease in pressure accompany- 
ing the reduction in Ccfl/V'* results in the disappearance of 
the pressure maxima and minima essential for two-phase 
equilibria. This is why four of the DLAs in Figures 12a 
and 12b are missing from Figure 12c. Clearly, more realis- 
tic models might include a range in Ccr/'4'* ratios rather 
than assigning the same value to each DLA. 

The point of this exercise is to show that the C II/C 
I ratio depends sensitively on the X-ray and cosmic ray 
ionization rates, both of which are uncertain. For these 
reasons, it is premature to use the C II/C I ratio to rule 
out CNM gas in DLAs. On the other hand, the observed 
C II/C I ratios are useful for placing upper limits on the 

Cci?/'0* ratio. We find that Ccfl/V"* in our model CNM 
gas cannot be more than 2 times the value in equation (9) 
in Paper I. Otherwise more than one third of the points in 
the Figure 12 would lie above the diagonal line. 

5.2. Equilibrium Gas Pressures in DLAs 

The second argument against CNM gas in high-z DLAs 
was made by Norman & Spaans (1997). They concluded 
that high-z DLAs instead consisted of pure WNM gas with 
pressures P < Pmin- Computing two-phase equilibria in 
the context of CDM models for galaxy formation, they 
found that at z > 1.5 the gas equilibrium pressure, P^q, 
exceeded the hydrostatic pressure, Phydro, in the mid- 
planes of protogalactic disks embedded in dark-matter ha- 
los. By contrast Peg was predicted to be less than Phydro 
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Fig. 13. — Plot showing equilibrium pressure, (i.e., 
{PminPmax)^^^), versus z and ncNM- Pressures computed for 
standard CNM model with "Gal" dust and minimal depletion, (a) 
Plot of P/k versus z. Two upper limits at logio(P/fc) < 2.2 K 
cm""^ correspond to DLAs along sightlines toward Q1108— 07 and 
Q2237— 06, which probably encounter only WNM gas. In that case 
the pressures would be higher than shown, (b) Plot of P/k versus 
ncNM- 



at z < 1.5. They concluded that the hydrostatic pressure 
available in model protogalaxies is insufficient to generate 
two-phase media at z > 1.5. The ratio Phydro/ Peg de- 
creases with increasing redshift because Peg increases due 
to the sharp decline of metallicity with redshift predicted 
by Norman & Spaans (1997; Spaans 2002). 

We tested these predictions using measurements of Ic to 
infer Peg. The results shown in Figure 13 were obtained 
with the CNM, "Gal" dust, minimal depletion model dis- 
cussed above. In Figure 13a we plot Peg versus z. The av- 
erage of the pressures corresponding to the positive detec- 
tions is Peg=(l-33±0.81)xl03 K cm'^. This is much lower 
than the gas pressures predicted by Norman & Spaans 
(1997), which exceed ~ 10^ K cm~'^ at 2 > 1.5 (see their 
Figure 3), but in agreement with Phydro predicted for typ- 
ical CDM models. According to the CDM models of Mo 
et al. (1998), 2/3 of DLAs detected in absorption at z 
« 2.5 should arise in halos with circular velocities, Vc > 
100 km s^^ . Using the Norman & Spaans (1997) formal- 
ism we find the corresponding hydrostatic pressures to ex- 
ceed 2x10'^ K cm~^. Given the uncertainties, we conclude 
that the hydrostatic pressures available in CDM models 
are sufficient to generate the type of two-phase equilibria 
we infer for DLAs. Coincidentally, these pressures could 
also arise in high-z DLAs if they resemble nearby mas- 
sive galaxies such as the Milky Way where the gas pres- 
sures at the solar circle are « 3x10^ K cm~^ (Wolfire et 
al. 2002). The analysis of Mo et al. (1998) and Nor- 
man & Spaans (1997) shows Phydro oc [V^7J(z)]^, where 
the Hubble parameter H{z) is an increasing function of 
redshift. Therefore low-mass objects at high redshift can 
in principle generate as much hydrostatic pressure as high- 
mass objects with low redshifts. As a result, pressure 
does not discriminate between DLA models based on CDM 
from the null hypothesis in which DLAs are the unevolved 
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disks of current normal galaxies (e.g. PW97). Nor docs 
the absence of redshift evolution evident in Figure 13a. 
This is an obvious prediction of the null hypothesis. In 
the case of CDM we note that the non-linear mass scale, 
Adc oc (l+z)^(*'/("+'^^\ where n is the power-spectrum 
index (Padmanabhan 1993), and the virial velocity, Vc 
oc [H{z)Mc{z)]^^^, where the Hubble parameter, H{z) oc 
(1 + z)^/^ for the Einstcin-dcSitter cosmology. Since n = 
— 2 at galactic scales, Phydro is independent of redshift. 
As a result, the redshift dependence of Phydro is negligible 
in both models. 

In Figure 13b we plot Peq versus CNM density, ticnm- 
The densities corresponding to positive detections (red 
data points) range between 2 and 100 cm~'^ with an av- 
erage < ncNM > = 16 cm~^. These resemble densities 
inferred for the CNM clouds in the Galaxy (W95). Be- 
cause the H I column density of a typical DLA velocity 
component is approximately 1x10^° cm~^ (PW97), the 
volume densities imply physical dimensions on the order 
of a few pc. This conclusion differs from arguments that 
DLA clouds causing 21 cm absorption uniformly cover the 
cores of compact radio sources, which typically have lin- 
ear sizes of ~ 100 to 400 pc (Briggs & Wolfe 1983), and 
has important implications for interpreting measurements 
of 21 cm absorption in DLAs (see § 5.3). The linear cor- 
relation evident in Figure 13b is tentative, because the 
Kendall r test for positive detections indicates r = 0.45 
and the probability for the null hypothesis of no correla- 
tion, pKendaii = 0.052. We find the average temperature 
< T >=190±130 K for this model. The temperatures are 
higher than diffuse clouds in the ISM due to the combina- 
tion of low densities and low metallicities. 

5.3. High Spin Temperatures at Large Redshifts 

The third argument against CNM gas at high redshifts 
stems from the high spin temperatures, Tg, deduced from 
21 cm absorption observations of DLAs. Whereas spin 
temperatures in nearby spirals are less than 300 K (Dickey 
& Lockman 1990), the spin temperatures in most DLAs 
exceed 500 K (Chengalur & Kanekar 2001). The discrep- 
ancy is greatest at 2 > 3 where several DLAs exhibit Ts 
> 2000 K. Because the kinetic temperature, T, equals the 
spin temperature in most scenarios, the high values of T, 
have been interpreted as indicators of gas in the WNM 
rather than CNM phase (Carilli et al. 1996; Kanekar & 
Chengalur 2001; Chengalur & Kanekar 2002). This poses 
an interesting dilemma, as C II* absorption, which must 
arise in the CNM in most DLAs, is detected in two DLAs 
with high inferred spin temperatures. 

How do we reconcile the spin temperature limit, Ts > 
4000 K (Ellison et al. 2001; Kanekar & Chengalur 2001), 
with the w 100 K temperature predicted for the CNM in 
the z = 3.387 DLA toward Q0201-I-11? To answer this 
question we consider the sizes of the CNM clouds rela- 
tive to the background radio source in Q0201+11. Follow- 
ing the discussion in § 5.2 we find that the linear dimen- 
sions of CNM clouds are less than 10 pc, since ucnm ~ 6 
cm~^ and we have assumed half the total H I column den- 
sity of 1.8x10^^ cm~^ to be CNM gas equally distributed 
among 5 or more velocity components. To estimate the lin- 
ear dimension subtended by the radio source at the DLA 
we note that VLBI observations at u = 1.6 GHz show 



Q0201-I-11 to subtend a solid angle less than wi. 6=2. 5x5.0 
mas^ (Hodges et al. 1984). The solid angle of the source 
must be larger than Wi.e at ly = 323.7 MHz, the frequency 
of redshifted 21 cm absorption, otherwise the brightness 
temperature of this radiation in the rest frame of the z 
= 3.61 QSO, Tb(1492 MHz)=1.7x lO^^K. This exceeds the 
lO^^K Compton limit restricting the brightness tempera- 
ture of sources such as Q0201-I-11, which belongs to the 
class of peak spectrum objects not exhibiting relativistic 
beaming (Phillips & Mutcl 1982). To be consistent with 
the Compton limit the source must subtend an effective 
linear diameter exceeding 40 pc at the DLA (where we 
assume f2M=0.3, nA=0.7, and h=0.7). As a result the 
typical CNM cloud in this DLA covers a small fraction 
of the background radio source. This contrasts with the 
conclusions of deBruyn et al. (1996) who did not consider 
the Compton limit. 

It is now possible to understand the discrepancy between 
Tg and the temperature of the CNM. In the optically thin 
limit the apparent 21 cm optical depth t(21)=/(7jvmX 
TcwM (21) where fcNM is the area covering factor of CNM 
gas, T(7iVAf(21) is the 21 cm optical depth of a CNM cloud, 
r(21)=ln[S'c/5„], and 5*^, and Sc are the 373.3 MHz flux 
densities at Doppler velocity v and in the continuum re- 
spectively. Chengalur & Kanekar (2000) placed a 1-a up- 
per limit of 0.011 on r(21). We shall be more conserva- 
tive and assume a 95 % confidence upper hmit of r(21) 
< 0.022, which reduces the lower limit on T, to 2000 K. 
Because half the sample DLAs show positive evidence for 
C II* absorption, while in 3 out of 25 DLAs there is prob- 
able evidence for WNM gas (see Paper I), we find 0.5 < 
fcNM < 0.88. We shall assume Jcnm = 0.67 as this is 
in accord with the observed relative occurrence of mul- 
tiple C II* velocity components in DLAs. As a result 
we find tcjvm(21) < 0.033. To compute NcNMiH I), 
the H I column density of the CNM cloud, we must esti- 
mate cr, the Gaussian velocity dispersion of the gas, since 
NcNMOil)= l.82xlO^^XTcNMi'2l)TcNMV2Tra. Assum- 
ing a = 9 km s~^ , the velocity dispersion corresponding 
to the velocity width to which the upper limits on t(21) 
apply, andTcNM « 100 K, we find A^cArM(HI) < 1.4x10^° 
cm~^. This is consistent with a total CNM H I column 
toward the optical continuum of 9x 10^° cm~^ distributed 
among 5 or more components, and a typical CNM HI col- 
umn of < 1.4x10^° cm~^ covering the radio source. In 
other words, the inferred spin temperature for the DLA 
toward Q0201-I-11 is high not because the temperature of 
the absorbing gas is high, but rather because the H I col- 
umn densities of the foreground CNM clouds are low. In 
fact there is evidence for a decrease in A^(H I) with in- 
creasing redshift. Pcroux et al. (2002) find the mean 
H I column density of DLAs, < N{m) > = 3.9x10^° 
cm-2 at z> 3.5, while < N{m) > = 8.2x10^° cm'^ at z 
= 2.4—3.5. Because these correspond to the total A^(H I) 
rather than the CNM portion of the gas, it is plausible to 
assume the mean H I column densities of CNM gas to be 
2x10^0 cm-2 at the DLA redshift, z = 3.38. Therefore, 
it is not improbable for iV(H 1)cnm to equal (1— 2)xl0^° 
cm~^ in this DLA. If this interpretation is correct, the high 
H I column density (iV(H I) = 1.8x10^1 cm'^) and 
C II* detected in absorption toward Q0201-I-11 indicate 
that a larger than average number of CNM clouds are lined 
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up toward the optical continuum source of this QSO. This 
is supported by the excessive number of low-ion clouds 
spread across a velocity interval of ~ 270 km s~^ (Ellison 
et al. 2001); by comparison the median velocity interval 
of DLAs is w 100 km s^^ (PW97). Note, Chengalur & 
Kanekar (2000) would not have detected this configura- 
tion, since it would have an apparent optical depth, t(21) 
= (10pc/40pc)^ x5xTcArM(21) < 0.01, which is less than 
the 2-(T upper limit of 0.022. 

Ellison et al. (2001) carried out an imaging study of the 
field surrounding Q0201-I-11 that provides a test of the 
CNM hypothesis. If C II* absorption in this DLA arises 
in WNM gas, the FUV intensities (at A = 1500 A), J^, 
would be higher than indicated by our CNM model. We 
find (J,y)M/AfM=3.8x 10^-^* ergs cm~^ s~^ Hz^^ sr^\ while 
("^i/)cJVM=4.3xlO~^^ ergs cm~^ s~^ Hz~^ sr~^ for the case 
of "Gal" dust and maximal depletion. Based on photomet- 
ric redshifts, the leading candidate for the galaxy respon- 
sible for DLA absorption is object G2, the high surface- 
brightness region of an i? « 25.3 galaxy separated by A6' 
= 2.9 arcsec from the QSO. If the FUV radiation emitted 
by G2 heats WNM gas responsible for C II* absorption, 
according to the bulge hypothesis G2 would detected in 
the i? band with flux density, S^a=4:TTJ„{A0)'^ / (I+zdla)^ , 
where v = {1 + znLA)^'o- In this case the predicted AB 
magnitude, i? = 21.3. If instead this WNM gas is heated 
by a uniform disk of sources centered on G2 and extend- 
ing across the QSO sightline, the surface brightness of the 
disk would be /Xii=26.4 mag arcsec"^. Both WNM sce- 
narios are ruled out by the Ellison et al. (2001) images. 
In the bulge scenario, G2 would be 4 magnitudes brighter 
than observed, while in the uniform disk scenario the sur- 
face brightness of the disk would be detectable at the 4-a 
level in regions on the side of G2 away from the QSO. On 
the other hand if G2 heats CNM gas responsible for 
C II* absorption, then R = 23.6 in the bulge scenario, 
while /i/j=28.7 mag arcsec^^ in the uniform disk scenario. 
Because of the large uncertainties in the R magnitude of 
G2, the bulge scenario is only marginally consistent with 
the data, while the uniform disk scenario is definitely con- 
sistent with the data. Therefore, the gas producing 
C ir absorption must be a CNM if it is heated by FUV ra- 
diation emitted by sources associated with G2. Of course, 
G2 may be incorrectly identified and the galaxy associated 
with the DLA could be located within the PSF of the QSO 
(Ellison et al. 2001), in which case C II* absorption could 
arise in either phase. This emphasizes the importance of 
obtaining a spectroscopic redshift for G2. 

More recently, Kanekar & Chengalur (2003) deduced a 
95 % confidence lower hmit of > 1.4x10'' K for the 
^;=3.062 DLA toward Q0336— 01. By contrast we infer an 
equilibrium temperature of « 100 K from the presence of 
C II* absorption in this DLA (see Table 1 in Paper I). 
In principle, we can explain the high inferred value of Ts 
with low values of iV(II 1)cnm in this object also. How- 
ever, much lower CNM filling factors are required due to 
the higher value of T^. We think it more likely that the 
high Ts inferred for this DLA is a byproduct of a complex 
radio-source structure. Though this source is observed to 
be compact at 5 GHz, it is likely to be extended at the 
much lower frequency predicted for redshifted 21 cm ab- 
sorption. In fact, it is possible that the radio emission 



is concentrated in two equally bright components, each 
of which is symmetrically displaced more than 10 pc from 
the optical continuum (see Phillips & Mutel 1982). In that 
case the gas detected in optical absorption need not inter- 
cept radio- frequency radiation along the line of sight. In 
any case it is difficult to understand how the lower limit on 
Tg can correspond to the temperature of neutral gas in any 
DLA, since 50 % of H is collisionally ionized at 1.5x10* 
K. 

6. SUMMARY AND CONCLUDING REMARKS 

This paper expands on a new technique for measuring 
SFRs in DLAs developed in Paper I. Namely, in Paper 
I we showed how detections of C II* 1335.7 absorption 
could be used to infer the SFR per unit physical area, tp^, 
in DLAs. We showed that a two-phase neutral medium, 
in which CNM gas is in pressure equilibrium with a WNM 
gas, is a natural byproduct of thermal equilibrium. We 
also found that C II* absorption lines detected in DLAs 
could arise in either phase, but that significantly higher 
values of ip^ were required if this absorption arose in the 
WNM. In this paper wc use cosmological constraints to 
show that while the line-of-sight likely encounters both 
phases, C II* absorption arises mainly in the CNM. Our 
specific conclusions are as follows: 

(1) We compute the SFR per unit comoving volume in 
DLAs, pif{z), by combining the mean SFR per unit area, 
< ip,t{z) >, with the observed number of DLAs per unit 
absorption distance. We obtain statistically significant re- 
sults for two redshift bins centered at z = 2.15 and z = 
3.70; these are the first quantitative measurements of p* 
at « 2. The results show (a) p*{z) for the WNM model 
is at least 10 times higher than for the CNM model, (b) 
p*{z) for the CNM model is in approximate agreement 
with independent determinations from luminosities mea- 
sured for flux-limited samples of galaxies, and (c) no evi- 
dence for evolution in the redshift interval, z = [1.6, 4.5]. 
We also compute the bolometric background intensities, 
Iebl, generated by the SFR histories, p*{z). In every case 
the WNM models predict Iebl above the observed 95 % 
confidence upper limits. By contrast, the Iebl predicted 
for the CNM models are consistent with this limit. As a re- 
sult, models in which C II* absorption arises in WNM gas 
are ruled out. Finally, we develop a "consensus" model, 
which accounts for the systematic errors arising from vari- 
ous model uncertainties. We also find < tp^:{z) > for DLAs 
appears to decrease significantly with decreasing redshift 
at ^; < 1.6 if DLAs evolve into ordinary galaxies. To com- 
pute < ^*(2:) > we average V* over Rhi rather than the 
de'Vaucouleurs radius. 

(2) Wc consider several consequences of our work. First, 
we evaluate the mass of stars and metals produced by the 
star formation histories we derived. Though the mass in 
stars produced by z = is consistent with the masses of 
current stellar populations in current galaxies, the mass of 
metals produced by z = 2.5 is significantly higher than is 
observed for DLAs at that redshift. Various solutions to 
this dilemma come to mind, including ejection of metal- 
enriched gas from DLA, making DLAs a transitory phase 
of galaxy evolution. But this results in an IGM metallic- 
ity, [M/H] — —1.2, which is a factor of 100 larger than 
observed in the Lya forest. Rather we favor a scenario 
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in which star formation and metal production occur in a 
centrally-focated bulge region displaced from the DLA gas. 
Second, we evaluate the bulge hypothesis and find that to 
within 10 % accuracy the predicted p^,{z) agrees with that 
predicted by the uniform-disk model used to derive the re- 
sults discussed above. From the empirical upper limit on 
Iebl, we find that the fraction of PUV radiation attenu- 
ated by optically thick dust in the bulge cannot exceed 0.7. 
Third, we search for evidence of connections between stars 
and gas. We find no evidence for a Kennicutt (1998) corre- 
lation between tp^ and iV(H I). This is naturally explained 
by the bulge hypothesis in which SFRs in the bulge region 
are unrelated to the H I column density in regions giving 
rise to DLA absorption. The lack of a Kennicutt relation 
is also consistent with the model in which star formation 
occurs in the same region creating DLA absorption. This 
is because i/'* is the SFR averaged over the entire DLA 
(5 kpc or more in most models), while N{li I) is sampled 
over a transverse dimension corresponding to the linear di- 
ameter of the QSO, i.e., ~ 1 pc. On the other hand when 
these quantities are averaged over the DLA sample, they 
are consistent with the Kennicutt relation, indicating that 
it may be present in a statistical sense. We also look for 
correlations between and metallicity, and and kine- 
matic velocity width. Marginal evidence for correlations 
were found in the cases of metallicity and low-ion veloc- 
ity width. Confirmation would indicate that both low-ion 
line width and metallicity are global parameters, which are 
determined by quantities such as dark-matter mass. The 
reasons for null correlations in the cases of and high-ion 
velocity width are the same as for the Kennicutt relation. 

(3) We discussed tests of the ideas presented here. First, 
we present statistically significant evidence for a correla- 
tion between Ic and k. This is consistent with grain pho- 
toelectric heating by DLAs with a limited range of SFRs, 
since Ic oc Keip^. We show that this correlation is not nat- 
urally explained by alternative heating mechanisms such 
as cosmic-ray ionization and X-ray photoionization. Sec- 
ond, we consider the possibility that the ^^3/2 and 
fine-structure states in C+ are populated by CMB radia- 
tion rather than collisions. In that case Ic would not equal 
the; c;ooling rate as we have assumed, hut woiild instead 
reflect the local temperature of the CMB. We test this 
hypothesis by comparing the observed Ic with {Ic)cmb, 
the spontaneous energy emission rate per H atom when 
CMB radiation alone populates the fine-structure states. 
We conclude that the CMB alone cannot explain the ob- 
served level of C II* excitation in 28 out of 30 DLAs. This 
is because the observed Ic do not show the sharp increase 
with redshift predicted for (Ic)cmb and because in all 28 
cases, Ic >> {lc)cMB- The two exceptions are DLAs with 
upper limits on 1^ and in which {Ic)cmb > 0.81c- The im- 
plied cooling rates are very low, and may indicate passage 
of the QSO sightlines through a WNM subjected to SFRs 
less than or equal to the mean <tp^> deduced for most of 
our sample. In the third test, we use the ratio of C II* 
to resonance-line optical depths, TZ{v) = r„(C II*)/T„(Si 
II), to probe the multi-phase structure of the gas. Specif- 
ically, TZ{v) in the WNM should be {Ic)wnm / {Ic)cnm 
times TZ{v) in the CNM. Because {Ic)wnm I {Ic)cnm ~ 
0.08, TZ{v) should decrease significantly at velocities cor- 



responding to the WNM. Such variations may have been 
detected in one DLA. We also show how the WNM is easy 
to hide, and we discTiss possible evidence for an increase 
in TZ{u) at velocities where V'* in the CNM increases. 

(4) We discuss possible objections to our picture of DLAs. 
The most persistent of these is the proposal that DLAs do 
not contain CNM gas, but rather are comprised only of 
WNM gas. The first argument against the CNM is related 
to the large C II/C I ratios detected in DLAs, since this 
ratio is predicted to be lower in the CNM of the Milky Way 
ISM. However, this conclusion is sensitive to the ratio of 
X-ray and cosmic-ray heating rates to the SFR per unit 
area. In our models the X-ray heating rate is considerably 
lower than in the ISM. For our model parameters we are 
able to construct reasonable models predicting C II/C I 
observations. The next argument against the CNM is that 
the hydrostatic pressures available in low-mass galaxy pro- 
genitors predicted for CDM theories are significantly lower 
than the values of Pmin predicted for two-phase models. 
In that case, CDM models cannot generate sufficient pres- 
sure to create a two-phase medium. We show this model 
is incorrect because the high values of Pmin were based 
on an underestimate of the DLA metallicities at large red- 
shifts. The third argument against the CNM stems from 
the high spin temperatures deduced from the low 21 cm 
optical depth, t(21) < 0.02, and high N{Yi I) (= 2x10^1 
cm~^) of a high- redshift DLA. How does one reconcile the 
large inferred spin temperature (T, > 2000 K) with the 
low temperatures (T ~ 100 K) indicated by the detection 
of C II* absorption in this DLA? This answer is related 
to the geometry and H I column densities, A^(H 1)cnMi 
of the CNM clouds. Specifically, our models predict CNM 
clouds with sizes less than 10 pc, which is smaller than 
the diameter of the background radio source. In that case 
we find iV(H 1)cnm ^ 10^° cm'^ for typical CNM clouds 
covering the high-z background radio source. Therefore, 
the inferred spin temperature is high not because the tem- 
perature of the absorbing gas is high, but rather because 
A'^(H T)cNM is low. As a result, the large value of N(R I) 
toward the optical continuum source implies a larger than 
average number of CNM clouds along this line of sight, 
which is observed. 

While we have attempted to evaluate important sources 
of systematic error, uncertainties still remain in our anal- 
ysis. For example, we cannot rule out the possibility that 
the grain size distribution is radically different from the 
MRN model used in our calculations. However, we think 
this is unlikely, as the MRN distribution nicely describes 
grain properties in local group galaxies, which are linked 
to high-z DLAs by a similarity in the comoving density 
of stars in the galaxies to the comoving density of neutral 
gas in the DLAs. While the composition of the dust is 
also uncertain, sensitive searches for the 2175 A feature 
will greatly reduce this source of systematic error. Sec- 
ondly, though we have presented evidence for a two-phase 
medium in our DLA sample, in particular we described ev- 
idence for th(> presence of a WNM in 3 DLAs, the thermal 
state of the gas may differ from the equilibria envisaged 
in the models. Heiles (2001) has argued that warm clouds 
in the ISM have temperatures of « 2000 K rather than 
the 8000 K predicted for the WNM. However, the C II* 
absorption we observe is unlikely to arise in such gas be- 
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cause at r > 1000 K the ratio of Ic to the total cooHng 
(i.e., heating) rates is only slightly larger than predicted 
for the WNM. As a result, SFRs similar to those predicted 
for the WNM would be required and would likely exceed 
the observed upper limits on the bolomctric background 
radiation intensity. There may be other caveats we have 
not discussed, but they are not obvious to us at this time. 

Finally, we discuss the similarity between pii{z) for DLAs 
and LBGs. Does it mean they are the same objects? 
Schaye (2001) suggested they were, and that DLAs were 
created by the intersection between the line of sight to 
the QSO and gas outflow from foreground LBGs. He ar- 
gued that since LBGs brighter than i? = 27 have comoving 
densities nco=0.016ft,^ Mpc~"^, they could account for the 
incidence of DLAs provided their H I cross-sectional radius 
r=19h~^ kpc for a cosmology with (f2M,f^A)=(0.3,0.7). 
With our measurement of we can compute the SFR 
per DLA, M*, because p*=M*nco- From Table 1 we find 
M*=40 Mq yr-i for h=0.7. SFRs this high are ruled out 
by results from Ha imaging surveys of DLAs (Bunker et 
al. 2001; Kulkarni ct al. 2001). Furthermore, in a recent 
search for a cross-correlation between DLAs and LBGs, 
Adelbcrgcr ct al. (2002) found tentative evidence that 
DLAs with R<25 were more weakly clustered than LBGs 
with i? < 25. As a result, current evidence suggests that 
DLAs and LBGs brighter than R=27 arc drawn from sep- 
arate parent populations. But more evidence is needed 
before this can be accepted as a robust result. If this 
is confirmed, the similarity between the SFRs per unit co- 
moving volume of LBGs and DLAs would be a coincidence. 
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